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FOREWORD 


The  BESRL  Work  Unit,  “Computerized  Models  for  the  Simulation  of  Policies  and 
Operations  of  the  Personnel  Subsystem-SIMPO-l“.  was  conducted  by  the  Statistical 
Reaaarch  and  Analysis  Division.  The  task  constituted  the  initial  undertaking  of  an 
oparations  research  requirement  described  in  the  Army  Master  Study  Program  under 
the  title,  “A  Simulation  Model  of  Personnel  Operations  (SIMPO)“  and  is  Project 
2Q066101M711.  “Army  Operations  and  Intelligence  Analysis."  under  the  auspices  of 
the  Army  Study  Advisory  Committee.  Sub-Work  Units  include:  a)  Operational  Analysis 
of  Personnel  Subsystems:  b)  Cataloging  and  Integration  of  Existing  Manpower  Models; 
c)  Development  of  Measures  of  System  Effectiveness;  d)  Development  of  Modeling 
Techniques;  e)  Design  and  Programming  of  SIMPO-I;  f)  Application  and  Evaluation  of 
Computerized  Models;  and  g)  Problem  Oriented  Language  for  Management. 


The  present  Technical  Research  Report  deals  with  the  development  of  a  computer 
model  which  can  be  used  in  evaluating  many  different  military  manpower  systems.  The 
model,  the  Qeneral  Matrix  Manipulator  (GMM).  is  based  on  the  movement  of  elements 
within  or  between  matrices  used  to  represent  the  different  personnel  categories  accord¬ 
ing  to  sipplied  rules.  The  two  dimensions  of  dte  matrices  represent  two  time  variables. 
Description  of  the  routines,  instructions  for  model  use,  and  a  simple  problem  are 
provided. 
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BRIEF 


Requirement: 

To  develop  a  generalized  mass-flow  model  which  can  be  used  to  simulate  many 
different  Army  manpower  systems  by  using  an  appropriately  coded  data  deck. 


Research  Product: 

The  GMM  is  a  modular  set  of  computerized  routines  which  move  elrrments  represent¬ 
ing  groups  of  persons  within  and  among  a  variable  number  of  matrices.  Two  time  vari¬ 
ables  can  be  depicted— for  example,  time  in  tour  and  time  in  service,  or  time  in  grade 
and  time  in  service.  Movement  is  controlled  by  priority-of-fill  cards  input  at  simulation 
time  and  may  be  changed  during  the  simulation  period. 


Utilization: 

The  GMM  has  been  used  in  developing  DISTRO,  a  SIMPO-I  model  required  by  PRIMAR 
II,  "Program  to  Improve  Management  of  Army  Resources,"  Monitor  Team  of  the  SIMPO 
Steering  Committee.  It  has  also  been  used  to  assist  in  the  evaluation  of  policies  on  the 
overseas  assignment  sequence  for  officer  personnel  and  on  the  reenlistment  of  members 
of  the  WAC  for  the  Office  of  Plans  and  Programs,  Office  of  Personnel  Operations. 
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The  U.  S.  Aony  Behavior  and  Systems  Research  Laboratory  Work  Unit, 
"Simulation  Models  of  Personnel  Operations  ^SIMPO),"  has  developed  a  series 
of  network  flow  models  designed  to  help  the  managers  of  Army  personnel  make 
better  estimates  of  accessions,  training,  promotion,  assignment,  and  dis¬ 
tribution  needs  and  capabilities  under  a  variety  of  policy  considerations. 
The  models  developed  are  of  two  general  types:  mass  flow  models  in  which 
groups  of  similar  persons  are  moved  through  a  series  of  states  in  a  de¬ 
terministic  manner,  and  entity  flow  models  with  which  individuals  and 
their  attributes  can  be  considered.  In  the  SIMPO  entity  models,  losses, 
gains,  transfers,  and  other  changes  are  stochastic  (uncertain),  the  out¬ 
come  of  each  simulated  event  depending  upon  the  value  of  a  random  number 
generated  by  the  computer.  In  the  entity  models,  a  number  of  different 
characteristics  can  be  treated  separately- -each  variable  can  be  permitted 
to  have  many  values  and  more  variables  can  be  considered.  In  bulk  flow 
models,  only  a  single  value  of  a  characteristic  can  be  represented  at  a 
node.  Additional  time-related  variables  are  monitored  by  using  an  array 
as  the  data  file  at  the  node. 

The  first  mass  flow  model,  DYNAMOD,  was  quite  general  in  concept. 

There  was  to  be,  at  each  node  represented  in  the  model,,  a  vector  of  num¬ 
bers  with  position  in  the  vector  representing  time  at  the  node.  During 
each  simulated  time  period,  the  subsets  of  individuals  then  completing 
an  assignment  were  transferred,  losses  were  taken,  the  system  vectors 
were  updated,  requirements  were  calculated,  and  new  assignments  were  made. 
Four  system  variations  were  modeled  using  this  concept,  producing  four 
separate  computerized  models  (1). 

New  problems  within  the  Army  rotation  system  became  so  urgent  with 
the  development  and  subsequent  replacement  of  large  forces  for  Vietnam 
that  the  DYNAMOD  concept  was  abandoned.  Additional  system  complexity 
was  modeled  with  matrix  nodes  and  provision  for  greater  flexibility  in 
the  selection  of  assignment  priorities.  Plans  were  made  to  develop  the 
SIMPO-I  "Grand"  Model.  The  resulting  model,  now  called  the  General 
Matrix  Manipulator  (®!M),  has  been  especially  designed  to  stress  flexi¬ 
bility.  It  is  particularly  useful  in  evaluating  the  one-time  problem  as 
contrasted  with  problems  that  must  be  periodically  reexamined  and  that 
can  be  handled  more  efficiently  by  models  specially  developed  to  simu¬ 
late  the  appropriate  system.  Examples  of  special  purpose  SIMPO-I  models 
are  DYROM  II  (2)  and  ACGMOD  (^).  Another  special  purpose  SIMPO-I  mass 
flow  model  has  received  wide  use  with  different  rotation,  training,  and 
reenlistment  problems.  The  Career- Noncareer  Model  (4)  is  a  rotation 
model  with  the  capability  of  evaluating  transition  from  first-term  to 
career  status.  The  three  special  purpose  models  are  more  efficient  for 
the  system  they  represent  than  the  GMM  would  be.  The  value  of  the  GMM 
lies  in  its  usefulness  with  diverse  systems  since  it  permits  the  analyst 
to  make  an  early  response  to  one-time  questions  on  the  effects  of  a 
policy  change. 


THb  SMM  MODEL 


This  generalized  mass  flow  model  was  developed  aroimd  the  concept 
of  maneuvering  elements  of  square  tabular  data  displays,  corresponding 
to  personnel  categories,  according  to  rules  input  at  the  time  of  model 
execution.  Number  of  tables  and  Interactions  between  the  tables  were 
to  be  flexible.  Possible  Interactions  were  to  be  from  a  limited  set  of 
actions:  all  or  a  portion  of  an  element,  or  several  elements,  of  a 
table  could  be  moved  to  the  corresponding  position  In  another  table,  or 
out  of  the  system  of  tables;  the  movement  could  also  be  a  shift  to  the 
first  column  and/or  row  of  another  table.  The  movements  were  to  be 
made  under  constraint  of  an  upper  limit  on  the  total  number  allowed  In  a 
table  and  a  lower  limit  on  the  time  of  stay  In  the  table.  Updating  the 
table  from  one  simulated  time  interval  to  the  next  was  to  Involve  shift¬ 
ing  each  element  of  the  table  to  the  next  row  below  and  the  next  column 
to  the  right.  System  renewal  elements  were  to  be  directed  to  the 
selected  table(s)  at  time  of  model  execution. 

It  was  thought  that  problems  In  the  functional  areas  of  accessions, 
assignment,  and  reassignment  Including  rotation,  distribution,  and  pro¬ 
motion  would  be  especially  amenable  to  evaluation  with  such  a  model. 

Since  the  actual  model  of  the  system  of  Interest  would  be  dependent  upon 
the  flow  rules  supplied  by  the  analyst/user  of  the  programs,  great  flexi¬ 
bility  for  modeling  different  flow  systems  would  be  offered. 

The  modeling  capabilities  actually  Included  In  the  present  system 
of  computer  programs  are  these: 

1.  Monitoring  of  two  different  time  measurements. 

2.  Variable  number  of  main  categories  of  personnel,  l.e.,  groups 
of  tabular  displays. 

3-  Variable  number  of  sub-categories  of  personnel  — separate 
tabular  displays. 

4.  Optl  .lal  rules  for  transferring  all  or  a  portion  of  elements  In 
one  table  to  another  table. 

Variable  sized  tables--up  to  48  rows  by  48  columns. 

The  first  of  these  capabilities,  monitoring  two  time  measurements, 
is  accomplished  by  calling  the  two  dimensions  of  the  square  table  (array) 
time  dimensions  and  by  updating  all  elements  In  the  array  each  time  the 
system  Is  advanced  from  one  simulation  period  to  the  next.  The  partic¬ 
ular  time  measure  represented  Is  optional  but  must  be  consistent  through 
all  parts  of  the  simulated  system.  Thus,  it  is  possible  to  simulate 
months,  quarters,  or  years  so  long  as  all  rates  supplied  are  consistent 
with  tha  same  time  period,  and  each  row  and  column  in  the  table  repre¬ 
sents  a  time  step  of  the  same  order.  Some  time  combinations  relevsnt  to 
Army  personnel  system  problems  are  tiroe-in-tour  and  time- in- system,  tirae- 
in-grade  and  time- in- system,  and  tiroe-ln-tour  and  tlme-in-enlistment. 
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The  next  tvo  capabilities  are  related  to  the  ntxnoer  of  classej  or 
groups  of  people  required  to  model  a  given  system.  An  example  of  main 
categories  vould  be  short  tour,  long  tour,  and  OONUS.  Sub-categories 
simulated  might  be  first,  second,  and  third  or  subsequent  short  tours 
within  the  short  tour  main  category;  and  before  short  tour,  after  one 
short  tour,  and  after  two  or  more  short  tours  within  the  other  two  main 
categories.  Another  possibility  might  involve  tour  areas  as  main  cate¬ 
gories  and  grades  as  sub- categories.  The  model  provides  the  option  of 
upper  limits  on  the  total  number  in  each  separate  sub-category.  The 
model  uses  a  symmetric  arrangement  of  sub-categories  within  main  cate¬ 
gories;  that  is,  if  one  main  category  requires  9  sub-catcgorles,  all 
main  categories  must  have  9.  (It  is  possible  to  have  dummy  sub-cate¬ 
gories.)  Since  the  model  is  disk  based,  the  number  of  separate  cate¬ 
gories  is  limited  by  disk  storage  area  rather  than  by  computer  main 
memory. 

Capability  4,  optional  rules  of  transfer,  is  the  most  important 
feature  of  the  model.  Flow  between  the  elements  of  any  two  tables  is 
possible.  The  order  in  which  the  model  considers  the  flows  can  be 
varied  in  any  manner  permitted  by  the  set  of  rules  used.  Since  each 
problem  modeled  requires  its  own  set  of  rules,  nv.ny  different  flow 
patterns  can  be  described,  making  it  possible  '.o  model  quite  different 
systems . 

The  fifth  capability,  variability  in  table  size,  makes  it  possible 
to  consider  different  limits  on  tour  lengths,  or  different  limits  on 
time  in  grade.  In  applying  the  size  variation,  a  full-sized  table  is 
carried  along  in  the  computer  by  the  model  but  a  limited  area  of  it  is 
used.  Maximum  table  size  is  48  rows  x  48  columns.  The  last  row  and 
column  are  used  tc  accumulate  elements  advanced  into  them.  Suppose  tour 
length  is  only  12  months.  The  model  would  carry  along  the  48  x  48  area 
but  would  function  within  the  12  +  1,  or  1^,  rows  and  48  columns  re¬ 
quired  by  a  12-month  tour.  A  rule  of  movement  would  need  to  be  supplied 
to  keep  the  l^th  row  empty  if  12  months  is  the  absolute  upper  limit  of 
time  allowed  in  the  tour. 


TECHNICAL  DISCUSSION  OF  MATRIX  MODELING 

The  GMM  may  be  thought  of  as  a  collection  of  matrix  nodes  connected 
by  a  flow  network.  Development  of  routines  that  could  be  used  flexibly 
to  maneuver  elements  or  combinations  of  elements  from  one  matrix  to 
another  led  to  the  examination  of  certain  classes  of  problems.  These  are: 

1.  Number  and  relationship  of  the  matrices 

2.  Transition  from  one  matrix  to  another 
Search  patterns 

4.  Updating 
Matrix  size 

6.  Capacitance 
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In  the  following  paragraphs,  these  problems  are  discussed  and  dis¬ 
tinctions  are  drawn  between  the  capabilities  of  the  SIMPO-I  GMM  and  that 
of  a  completely  general  matrix  manipulator. 


Number  and  Relationship  or  Matrices 

In  a  personnel  subsystem,  there  Is  usually  a  constraint  on  total 
number  of  members.  This  limit  may  apply  to  the  lowest  subset  of 
Indlvldual8--a8  In  the  element  of  a  matrjx--or  It  may  apply  to  the  sum 
of  several  categories,  the  total  number  within  a  node  or  across  several 
nudes.  In  a  matrix  manipulator.  It  should  be  possible  to  structure  the 
.sysrem  under  either  type  of  constraint.  Either  one  or  more  matrices  may 
oe  allowed  In  a  cluster,  and  either  the  Individual  matrix  or  the  cluster 
may  be  capacitated.  The  two  matrix  '.Umenslons  may  be  as  large  as  the 
particular  system  reaulres  but.  In  the  Interest  of  reducing  slmu'atlon 
time,  should  be  kept  as  small  as  possible.  Also,  the  matrices  m^y  be  as 
numerous  as  required  by  the  system  under  consideration,  limited  only  by 
disk  area  available  and  simulation  time  costs. 


Transition  from  one  Matrix  to  Another 

In  a  time-dependent  model,  transition  from  a  cell  of  one  matrix  to 
another  cell  of  a  second  matrix  can  be  accomplished  by  one  or  another  of 
several  flow  types .  Denoting  the  losing  matrix  node  as  C  and  the  gain- 
Ing  matrix  node  as  D,  the  flow  types  In  terms  of  the  element  In  D  which 
may  gain  from  an  element  In  C  are  shown  In  the  following  table: 


Flow  Type 
1 

2 

5 

4 


Element  losing 
In  C  at  time  t 


Element  gaining  In 
D  at  time  t  +  1 


'ij 


1+1,  J+1 


ij 


1,  j+1 


'ij 


1+1,  1 


'11 


The  nodes  in  a  GMM  can  be  thougat  of  as  being  located  on  a  grid  on 
which  locations  are  specified  by  tour  and  level  (grade,  rank,  or  skill). 
When  more  than  one  node  is  located  at  a  grid  location  (e.g.,  when  several 
MOS's  are  Interchangeable  In  one  tour  but  not  necessarily  In  other  tours). 


these  nodes  are  considered  to  be  subtours  of  the  same  tour.  The  flow 
upward  to  higher  levels  within  a  tour  constitutes  the  familiar  feeder 
pattern,  the  flow  across  tours  at  the  same  level  depicts  the  rotation 
phenomenon . 

Thus,  the  GMM  is  a  collection  of  nodes  having  one  kind  of  flow 
(Type  1)  within  a  command  aggregation  (the  feeder  pattern)  and  a  second 
kind  of  flow  (Type  2)  across  comnands.  These  two  kinds  of  flow  are  de¬ 
picted  in  Figure  1  in  connection  with  Type  A  nodes  (matrices  of  time  in 
command  vs  time  in  service).  The  two  kinds  of  node  envisaged  for  the  GMM 
are  shown  in  Figure  2.  The  second  type  of  node  (Type  b)  consists  of 
matrices  whose  rows  represent  time  in  command  and  whose  columns  represent 
time  in  grade  or  skill  level.  Only  one  type  of  node  will  usually  be 
present  in  a  single  flow  pattern,  although  a  third  dimension  could  be  de¬ 
picted  by  Increasing  the  number  of  matrices,  in  effect,  combining  the 
characteristics  of  A  and  B  nodes  in  the  same  model. 

Type  1  flows  involve  adding  the  amount  taken  from  a  feeder  node  to 
the  equivalent  cell  of  the  gaining  node.  In  Type  nodes,  which  are  con¬ 
sidered  in  the  discussion  that  follows,  neither  time- in- service  nor  time- 
in-command  changes  as  a  result  of  the  flow.  Prescribed  proportions  can 
be  taken  from  each  row  (representing  those  with  equal  time  in  service). 

Type  2  flow  is  Implemented  by  adding  K  persons  to  a  cell  in  the 
gaining  node  that  corresponds  to  the  first  column  (least  time  in  cottinand 
and  the  row  (time  in  service)  from  which  the  K  persons  were  removed  in  the 
contributing  node.  Each  cell  in  the  matrix  would  thsn  be  advanced  one 
row  and  one  column.  Nondeployability  constraints  relating  to  rotation 
status  are  imposed  on  this  flow.  In  a  distribution  model  having  both 
Type  1  and  Type  2  flows,  more  detailed  distribution  of  assets  to  claimants 
could  undoubtedly  be  made  after  the  Type  1  flows  and  before  the  Type  2 
flow.  The  rate  of  fill  (constraints  on  flow  into  nodes)  at  each  connand 
aggregation  could  then  be  recomputed  to  reflect  distribution  priorities. 

A  Type  3  flow  would  represent  a  flow  within  a  comnard  via  the 
feeder  pattern  among  Type  B  nodes.  These  nodes  have  time  in  node 
(grade  or  skill  level  for  an  MOS  or  branch)  as  the  vertical  dimension  of 
the  matrix  and  time  in  command  (tour)  as  the  horizontal  dimension.  The 
Type  3  flows  will  have  the  effect  of  resetting  time  in  node  to  zero  and 
retaining  the  column  location  (time  in  coninand)  of  the  receiving  cell. 

This  flow  could  be  used  to  depict  an  officer  promotion  policy. 

Type  4  flows  are  included  here  for  the  sake  of  completeness.  The 
Type  4  flows  would  place  all  input  into  a  node  at  the  upper  left-hand 
cell  (i.e.,  with  time  in  grade  and  time  in  command  both  set  at  zero). 

For  a  Type  B  node,  this  flow  would  represent  promotion  to  fill  a  vacancy 
located  in  a  different  command.  See  Figure  3* 
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Figure  3.  Node  Type  B  flows  in  the  GMM 


Search  Patterns 

Flows  between  model  nodes  are  under  the  control  of  priorlty-o£-£lll 
rules  provided  by  the  analyst.  The  order  In  which  nodes  are  reached  in 
extracting  people  to  fill  other  nodes  is  determined  by  the  sequence  of 
the  rules.  Another  aspect  of  patterns  of  search  concerns  the  order  in 
which  elements  of  the  same  losing  matrix  are  entered  and  removed.  To 
simplify  development  of  rules  of  flow  in  the  manipulator,  it  is  desirable 
that  areas  of  the  losing  matrix  be  covered  automatically  by  specifying 
appropriate  transfer  limits  in  terms  of  rows  and  columns  (time)  beyond 
which  the  cells  of  a  particular  matrix  are  not  decremented  to  provide 
flow  from  the  node.  Direction  of  extraction  from  the  transfer  area  is 
also  important.  In  one  Instance,  the  search--and  the  successive 
emptying- -begins  with  the  highest  numbered  row  (column)  and  works  down, 
while  in  the  other,  search  starts  with  the  lowest  numbered  row  and  works 
up.  Complete  flexibility  requires  Chat  the  model  user  be  given  the 
option  of  covering  specifically  any  of  the  shaded  areas  of  the  matrices 
in  Figure  4,  with  search  progressing  in  either  direction.  All  these 
search  areas  can  be  covered  by  specifying  upper  and  lower  row  and  column 
boundaries  (matrix  (n)  in  the  figure).  Similarly,  direction  of  search 
can  be  determined  by  specifying  increasing  or  decreasing  dimensionality 
for  each  matrix  dimension. 


As  is  often  the  case  when  plans  for  a  model  are  being  developed, 
the  concepts  of  search  patterns  discussed  above  were  not  immediately 
apparent  to  the  model  developer.  The  SIMPO-I  Matrix  Manipulator  makes 
possible  flow  Types  1,  2,  and  3  with  all  searches  starting  at  the  upper 
time  limit  and  progressing  toward  the  lower. 


Updating 

Bringing  time  up  one  unit  is  a  different  problem  in  the  mass  flow 
models  than  in  the  entity  models.  In  the  mass  flow  models,  position  in 
the  node  array  indicates  tire  spent  in  the  node,  node  cluster,  or  system. 
In  the  entity  models,  time  in  state  is  computable  from  state  entry  time 
and  present  simulation  time.  Moving  the  nth  element  of  the  mass  flow 
node  array  to  the  (n+l)th  position  poses  no  problem  until  the  array 
boundaries  are  encountered.  This  problem  can  be  handled  in  the  fixed 
length  node  by  making  certain  the  priority-of-fill  rules  move  the 
element  into  the  last  position  before  updating  takes  place — or  as  a 
first  step  of  the  updating  process.  The  node  representing  residual 
system  members  not  in  fixed-length  nodes--or  not  in  a  high  priority 
function- -must  either  have  safeguards  that  prevent  the  loss  of  the  last 
element  at  updating  or  the  node  array  must  be  so  large  that  the  last 
position  always  remains  empty.  In  the  GMM,  the  last  element  in  each 
row  or  column  is  cumulative,  and  the  time  represented  for  this  part  of 
the  matrix  is  "greater  than  t"  where  t  is  the  nominal  duration  of  the 
node  on  the  single  dimension. 


Uniform  or  Different-Sized  Matrices 

The  updating  function  •'s  related  to  another  basic  concept  for  a 
matrix  manipulator--that  of  size  of  array  for  each  node.  In  modeling  a 
personnel  system,  some  nodes  can  be  easily  classified  as  a  matrix,  othert 
are  rightly  a  vector,  still  others  an  element  (a  pool).  The  Phase  I  GMM 
uses  a  matrix  at  each  node.  The  matrices  are  square  and  all  are  the  samt 
size.  Within  the  square  matrix,  it  is  possible  to  use  a  second  limit 
(length  of  tour)  on  the  number  of  rows  in  which  an  element  other  than  a 
zero  may  appear.  Suppose  the  OONUS  tour  is  24  months  long  and  the  maxi¬ 
mum  matrix  size  is  48  by  48.  With  24  and  48  input  as  model  parameters, 
individuals  moving  past  24  months  in  a  CONUS  assignment  are  accumulated 
in  the  23th  row.  Updating  on  the  other  dimension  continues  during  each 
subsequent  period  until  the  48th  position  in  the  23th  row  is  reached,  at 
which  point  further  accumulation  occurs.  The  effect  is  that  the  modeler 
can  stipulate  a  23  by  48  matrix  at  a  particular  node. 

At  the  present  time,  it  seems  desirable  to  modify  the  existing  mode 
to  include  the  capability  of  handling  variable  dimensions  for  the  nodes. 
The  Phase  I  GMM  uses  a  combination  of  maximum  matrix  size  and  matrix 
number  sequence  to  reach  the  correct  position  on  the  disk  where  a  matrix 
node  is  written.  Changing  the  concept  of  the  model  from  a  symmetrical 
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arrangement  of  sets  and  subsets  of  square  matrices  to  a  nonsymmetrlcal 
one  composed  of  elements^  vectors,  and  matrices  as  required  by  the  sys- 
ten  being  simulated  may  mean  major  changes  in  all  computer  routines  or 
complete  restructuring  of  the  model.  While  the  SIMPO-II  GMM  will 
probably  use  variably  dimensioned  nodes,  it  is  well  to  remember  that 
the  SIMFO-I  GMM  uses  a  square  matrix  at  each  node  and  accesses  the  disk 
by  a  combination  of  dimension  and  matrix  sequence  number. 

In  orde'*  to  circumvent  meaningless  search  of  empty  cells  of  the 
losing  matrix,  the  SIMPO-I  GMM  keeps  a  record  of  the  extreme  matrix 
dimension  at  which  non-zero  elements  are  located.  Search  starts  at  the 
extreme  non-zero  element  and  works  back  toward  row  or  column  number  one. 


Matrix  Capacitance 

Jus*"  as  provision  must  be  made  to  fix  an  upper  limit  on  the  number 
of  persons  contained  in  a  given  matrix  node  (or  set  of  matrix  nodes),  a 
truly  genetal  model  must  also  consider  the  lower  limit.  In  the  present 
GMM,  the  lower  limit  if.  considered  only  in  the  special  purpose  DISTRO 
application,  where  distribution  of  available  personnel  to  Army  elements 
within  a  group  of  matrices  is  simulated  by  using  a  special  computer 
routine  with  results  of  a  GMM  simulation  period  as  input.  Later  develop¬ 
ments  (SIMPO-II)  in  the  GMM  will  involve  lower  as  well  as  upper  limits 
for  matrix  capacity. 


Model  Input  and  Output 

The  GMM  requires  the  starting  data  to  be  supplied  to  the  model  in 
the  same  configuration  in  which  it  is  to  be  used.  In  simulating  a  com¬ 
plicated  persormel  system,  data  preparation  is  a  tedious  clerical  task 
even  when  appropriate  information  is  available.  The  necessary  details 
are  often  not  av'ailable  in  a  form  which  allows  the  analyst  to  make  the 
necessary  suimnsries.  Since  the  Army  keeps  its  personnel  records  on 
tapes  and  inak«vs  extracts  available  for  operational  and  study  needs,  it 
is  feasible  to  provide  routines  which  allow  the  GMM  to  obtain  data  with¬ 
out  requiring  intermediate  hand  processing.  These  routines  are  planned 
as  a  part  of  the  SIMPO-II  effort. 

Planning  appropriate  model  output  for  a  general  matrix  manipulator 
is  more  difficult  than  for  a  specific  model.  System  status  for  a  speci¬ 
fied  time  period  is  perhaps  the  most  easily  provlded--and  the  most 
difficult  to  Interpret,  unless  the  simulated  system  is  a  very  simple  one. 

A  summary  based  on  flows  between  nodes  or  on  the  total  number  of  persons 
in  a  given  matrix  node  at  a  specified  period  is  probably  more  useful 
information.  The  Gl-IM  package  contains  routines  to  obtain  these  two  types 
of  aggregated  information.  It  is  also  possible  to  sum  the  elements  in  a 
given  area  of  one  or  more  matrices.  For  example,  it  would  be  possible  to 
obtain  the  sum  of  all  persons  who  have  been  in  short  tour  less  than  6 
months,  .f  time  in  tour  were  one  of  the  time  variables  und^tr  consideration. 
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This  capability  is  useful,  too,  when  the  first  month  of  an  assignment 
area  is  to  be  considered  transient  time  and  duty  personnel  are  obtained 
from  those  there  longer  than  a  month. 


SIMPOl  GMM  ROUTINES 

To  insute  flexibility  and  adaptability  to  special  problem  areas, 
the  GMM  is  conceptualized  as  a  network  of  matrix  routines  interwoven 
into  a  logical  unit.  As  new  specifications  arise,  this  network  of 
matrix  routines  can  be  expanded  to  include  additional  special  function 
routines.  The  only  major  change  necessary  within  the  GMM,  besides  the 
development  of  special  routines,  will  be  the  logical  integration  of  the 
new  routines  into  the  main  GMM  driver  program.  This  network  concept  has 
not  completely  materialized  in  the  SIMPO-I  version  of  the  GMM,  for  the 
current  main  GMM  program  accomplishes  some  functions  which  will  be  accom¬ 
plished  by  separate  routines  in  the  SIMPO-II  GMM. 

One  main  routine  serves  as  a  driver  program  for  the  entire  GMM;  it 
either  performs  a  specific  function  or  calls  the  appropriate  subroutine 
which,  in  turn,  performs  this  function.  As  a  result,  the  main  routine 
determines  the  logical  basis  of  the  simulation,  and  the  validity  of  any 
model  is  directly  related  co  the  logic  employed.  The  main  routine  driver 
and  its  Interface  with  the  individual  subroutines  is  described  in  the 
following  sections. 


Main  Oriver  Program,  MAINGMM 

MAINGMM,  the  driver  program,  is  a  logical  group  of  statements  which, 
using  an  iterative  procedure,  simulates  each  tl'ie  period.  These  state¬ 
ments  define  the  simulation  logic  by  the  order  and  manner  in  which  t!  t-y 
employ  the  separate  GMM  subroutines. 


MAINGMM  logic  can  be  sumnarized  in  a  series  of  steps,  each  repre¬ 
senting  a  special  function  or  event  within  the  simulation  process. 
These  steps  perform  the  following  functions: 


1. 

2. 

5. 

4. 

5. 

6. 

7. 

8. 

9. 


input  simulation  data 
update  the  system 
make  initial  transfers 
calculate  node  assets 
determine  node  requirements 
fill  node  requirements 
make  final  transfers 
input  additional  personnel 
increase  time  period 


These  steps  are  accomplished  within  MAINGMM  in  conjunction  with  the 
separate  subroutines.  In  the  SIMFO-I  GMM,  the  MAINGMM  driver  per¬ 
forms  steps  1,  2,  8,  and  9  without  the  aid  of  subroutines. 
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Figure  5.  Data  Deck  Setup  for  GMM 
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Additional  subroutines  will  take  over  these  five  functions  In  the 
SIMPO-II  phase  of  the  GMM.  Subroutines  TIMEUP,  SUMMARY,  FILUP, 

OUTSIDE,  MX,  X^B,  and  LOCK  accomplish  steps  3,  <ind  7.  Each  of 

these  nine  steps  and  the  Individual  subroutines  are  described  In  detail 
below  In  the  order  in  which  they  occur  In  the  Iterative  process. 

Step  1.  Input  simulation  data.  The  '*ata  Input  deck  (Figure  5)  can 
be  divided  Into  three  sections,  each  defined  according  to  Its  function 
within  the  simulation.  These  three  Input  sections,  outlined  In  Table  1 
and  described  In  Table  2,  determine  respectively:  the  status  of  the 
system,  the  directional  flow  of  personnel  within  the  system,  and  the 
special  monthly  system  requirements. 

Status  of  the  system.  The  primary  parameter  card  and  the  tcur  deck 
setup  describe  In  detail  the  system  that  will  be  simulated.  The  primary 
parameter  card  sets  the  boundary  conditions  within  which  the  system 
operates.  That  Is,  It  determines  the  number  of  nodes  and  node  clusters, 
the  number  of  flow  rules,  and  specific  program  options. 

The  tour  deck  setup  cards  describe  each  node.  Its  characteristics, 
and  personnel  within  that  node.  Tour  deck  card  A  data  serve  as  parameters 
for  each  node,  defining  the  node  In  relation  to  other  nodes  and  describing 
Its  length  and  loss  rate.  Tour  deck  card  B  data  describe  personnel  dis¬ 
tributed  within  the  node.  Currently,  an  accurate  distribution  of  personnel 
according  to  time  served  within  a  particular  node  Is  practically  Impossible 
to  obtain.  Therefore,  a  uniform  distribution  of  time  served  in  the  node  Is 
usually  assumed. 

Since  the  tour  deck  setup  card  B  data--the  personnel  dlstrlbution-- 
comprlses  the  bulk  of  the  GMM  data,  card  B  data  are  stored  on  an  on-line 
disk  to  save  computer  storage  space.  These  data  can  easily  be  retrieved 
as  needed.  Inputting  data  In  segments  onto  a  permanent  storage  disk  is 
definitely  more  practical  than  reading  all  the  data  from  cards  Into  the 
core  storage  area,  saving  computer  time  as  well  as  computer  storage  space. 

The  MAINGMM  has  an  option  to  read  the  card  B  data  from  either  a 
permanent  storage  disk  or  from  the  cards  directly.  If  the  data  have 
previously  been  stored  on  the  permanent  disk,  the  user  Inputs  only  ITT 
cards  A;  the  program  atuomatlcally  retrieves  the  card  B  data  from  the 
permanent  storage  disk.  Only  If  the  user  turns  on  sense  switch  3  wHl 
the  program  read  both  the  card  A  and  card  B  data  from  cards.  In  either 
case,  as  the  card  B  data  are  input  In  segments  to  the  storage  area 
within  core,  card  B  data  are  transferred  to  a  temporary  working  storage 
disk.  Since  sev’eral  computer  runs  are  usually  made  to  compare  output 
from  the  different  model  options,  the  permanent  disk  storage  prevents 
duplicative  reading  of  the  same  data  cards  onto  the  temporary  working 
disk.  Instead,  the  data  can  be  read  directly  from  the  permanent  storage 
disk  and  transferred  to  the  temporary  disk. 
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Table  1 


INPUT  TO  OIM 


No.  of 

Card  Description 

Cards 

Variable 

Columns 

Format 

Primary  Parameter 


System  Status  Parameters 
1 


(M  *  ((MAXLEN  -  l)/20 


ITT  Cards  A. 


Tour  Deck  Card  A 


ID 

1-5 

A5 

NTOUR 

4-7 

14 

ITT 

8-11 

14 

FIRST 

12-15 

14 

LAST 

16-19 

14 

IFILL 

20-23 

14 

NP 

24-27 

14 

MAXSUB 

28-31 

14 

CIOS 

32-35 

14 

JC 

56-39 

14 

IC 

40-43 

14 

MAXLEN 

44-47 

14 

LRT 

48-51 

14 

PDW 

52-57 

F6.4 

IDISTON 

58 

11 

I PUNCH 

59 

11 

(found  in  columns  8-11  of  Primary 
with  one  Card  A  followed  by  M  Cards  B 

runcated  *  MAXLEN) 

Sense  Switch  is 

• 

on.  Otherwise 

input 

TYPE 

1-10 

110 

SUB 

11-20 

110 

ACTUAL 

31-40 

110 

NUM 

41-50 

no 

LENGTH 

51 -60 

no 

OUT 

61-66 

F6.3 

Construct  M  (M  -  ((MAXI.EN  -  l)/20  »  l)  truncated  *  MAXLEN)  Cards  B: 
(MAXLEN  is  found  in  columns  44-47  of  the  Primary  Parameter  Card). 
This  means  the  matrix  data  are  input  row  by  row  with  MAXLEN  rows  and 
(MAXLEN  -  l)/20  +  1  truncated  cards  per  row. 
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Table  1  (Continued) 


Card  Description 

No.  of 
Cards 

Variable 

Columns 

Format 

Tour  Deck  Cards  B 

M 

SYST(1,1) 

1-4 

14 

SYST(1,2) 

• 

5-8 

14 

• 

• 

• 

SYST(1,20) 

77-80 

• 

• 

14 

SYST(1,21) 

• 

1-4 

14 

• 

• 

• 

SYST(1,40) 

77-80 

• 

• 

14 

SYST(1,41) 

• 

1-4 

14 

• 

• 

* 

SYST{l,M/dCLEN) 

• 

• 

14 

SYST(2,1) 

• 

1-4 

14 

• 

• 

« 

SYST(MAXLEN, 

• 

• 

14 

maxlen) 

Directional 

Flow  Parameters 

Transfer  Rule 

1 

NOTT 

1-4 

14 

Parameter  Card 

NOFIRST 

5-8 

14 

NOLAST 

9-12 

14 

Transfer  Rules 

NOTT 

OUTTT 

1-10 

no 

OUTST 

11-20 

no 

OUTTO 

21-50 

no 

OUTSO 

51-40 

no 

PERD 

41-50 

no 

PEROTO 

5I-6O 

no 

PCT 

61-66 

F6.5 

Prior ity-of- fill 

NP 

INTOUR 

1-8 

18 

Rules 

INSUB 

9-x6 

18 

OUTTOUR 

17-24 

18 

OUTSUB 

25-52 

18 

AFTER 

55-40 

18 
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Table  1  (Continued) 


Card  Description 

No.  of 
Cards 

Variable 

Columns 

Format 

PEROUT 

41-46 

F6.5 

PER 

47-52 

F6.5 

REP 

55-58 

F6.5 

lOD 

59-61 

15 

IGRADE 

62-64 

15 

ITYPE 

65-67 

13 

Last  Resort  Tour 

LRT 

LSTRSTO 

1-5 

15 

Rules 

LSTRSTT 

6-10 

15 

LSTRSTS 

11-15 

15 

Monthly  Requirements 

Monthly  requirements  consist  of  one  group  of  two  sets 

of  cards. 

one  set 

of  N  Quota  cards  and  one  set  of  M  Input  cards, 

for 

each  time 

!  period  to 

be  simulated. 

(n  -  ((ntour  - 

l)/10  +  1  truncated;  M  ■  ((CIOS 

-  1)/10  +  1) 

truncated) . 

Monthly  requiresments  consist  < 

of  one  group  of  two 

cards,  one  Quota  card 

and  one  Input 

card,  for  each  time  period  to  be 

simulated 

Quota  Card 

N 

NEEDS(l) 

1-8 

18 

NEEDS(2) 

• 

9-16 

18 

• 

• 

• 

NEEDS(  NTOUR) 

• 

18 

Input  Card 

M 

lOS(l) 

1-8 

18 

I0S(2) 

e 

9-16 

18 

• 

« 

• 

lOS(CIOS) 

• 

18 

16 


Table  2 


VARIABLES  INPUT  BY  USER  TO  GMM 


Variable 


Definition 


ID 

NTOUR 

MAXSUB 

ITT 

FIRST 

LAST 

IFILL 


NP 

MAXLEN 

CIOS 

PDU 


Simulation  Identification  (3  characters) 

Number  of  node  clusters 

Maximum  number  of  nodes  within  each  node  cluster 
(No.  of  subtours) 

Maximum  number  of  nodes  (NTOUR  *  MAXSUB) 

First  period  to  be  Simulated 

Last  period  to  be  simulated 

Parameter  governing  filling  of  requirements: 

If  IFII.L  leas  than  1 — have  NTOUR  requirements,  one  for 
each  node  cluster. 

If  IFILL  .GE.  l--have  ITT  requirements,  one  for  each 
node. 

Number  of  priority  of  fill  rules 

Maximum  length  of  system  in  time  periods 

Number  of  categories  input  from  outside  the  system 

Percentage  of  requirements  which  may  be  filled  from 
outside  the  system 


ITT  of  the 

TYPE(I)* 

SUB(I) 

ACTUAL(I) 

NUM(I) 


following: 

Node  cluster  designation  (Numbered  one  to  NTOUR) 
Node  category  designation  (Numbered  one  to  MAXSUB) 
Number  of  entities  in  node 
Number  of  nodes  within  the  node  cluster 


LENGTH(I) 

OUT(I) 

NOTT 

NOFIRST 

NOLAST 


Length  of  node  in  time  periods 

Percentage  of  entities  lost  from  node  at  end  of  12  time 
periods  (yearly  attrition  rate) 

Number  of  possible  horizontal  node  to  node  movements 
( transfers) 

Number  of  initial  transfers  or  horizontal  movements 
Number  of  final  transfers  or  horizontal  movements 
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Table  2  (Continued) 


Variable 


Definition 


NOTT  of  the  following: 

Node  cluster  to  which  horizontal  movement  iomade  (gaining  cluster ) 
Node  within  node  cluster  to  which  movement  is  made  (gaining  node ) 

Node  cluster  from  which  horizontal  movement  is  made  (losing  cluster) 
Node  within  node  cluster  from  which  movement  is  made  (losing  node) 


OUTTT( J) 
OUTST(J) 
OUTTO(J) 
OUTSO(J) 
PERD( j) 


Row  of  node  (particular  time  period)  from  which 
horizontal  movement  is  maae 

PERDTO(j)  Row  of  node  to  which  horizontal  movement  is  made 

PCT(j)  Percentage  of  entities  within  each  cell  of  losing  node  which 

are  available  for  horizontal  movement 


NP  of  the  following: 

INTOUR(K)'  Node  cluster  to  which  movement  is  made  (gaining  cluster) 

INSUB(K)  Node  within  node  cluster  to  which  movement  is  made  (gaining  node) 

OUTTOUR(K)  Node  cluster  from  which  movement  is  made  (losing  cluster) 

OUTSUB(K)  Node  within  node  cluster  from  which  movement  is  made  (losing  node) 

AFTER(K)  Row  (time  in  tour)  of  node  from  which  movement  is  to  be  made 

(see  explanatory  note  in  IOD(k)) 


PER(K) 


PEROUT(K) 


lOD(K) 


REP(K) 


Percentage  of  requirements  of  node  to  which  movement  is 
possible  which  may  be  made  from  the  node  cluster 
(OUTTOUR)  and  node  (OUTSUB)  specified 

Percentage  of  node  (OUTTOUR)  that  may  be  used  as  fill 
for  ( INTOUR )  node  to  which  movement  is  specified 
(Assumed  to  be  100^  if  blank  on  input) 

If  lOD  .EQ.  l--movement  is  only  from  time  period  (row) 
Indicated  (AFTER(K)) 

If  lOD  .NE.  l--moveraent  is  from  time  period  MAXLEN 
through  AFTER(K) 

i.e.  if  AFTER  .EQ.  12,  movement  takes  place  from  row  15  only 
if  lOD  .EQ.  1;  if  lOD  .NE.  1,  from  MAXLEN  to  12th  rows 

If  REP  .EQ.  l--may  fill  NEEDS  with  100^  of  those 
available  from  outside  the  system 

If  REP  .NE.  1--NEEDS  are  filled  with  PEW^  of  those 
available  from  outside  the  system 
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Table  2  (Continued) 


Variable  Definition 

ITYPE(K)  Determines  how  data  will  be  input  into  a  node 

If  ITYPE  .LE.  2 — input  will  be  into  row  1,  column  T. 

If  ITYPE  .EQ.  ^''input  will  be  into  row  1,  column  1. 

If  ITYPE  .GE.  4 — input  will  be  into  row  1,  column  1. 

LRT  Number  of  categories  of  personnel  available  as  input  to 

system  which  must  be  emptied  and  saved  in  the  event 
they  are  not  needed  by  algorithm 

LRT  of  the  following: 

LSTRSTT(L)*  Tour  to  which  LSTRSTO(L)  input  is  sent 

LSTRSTO(L)  Designation  of  cell  containing  input  to  system 

LSTRSTS(L)  Node  within  node  cluster  to  which  LSTRSTO(L)  input 
is  to  be  sent 


NEEDS(M)* 

NE(I)* 

lOS(N)' 


If  IFILL  .EQ.  0--Read  in  NTOUR  node  cluster 
requirements  for  the  period  being  simulated 

If  IFILL  .EQ.  1--Read  ITT  node  requirements  for 
period  being  simulated 

Number  input  from  outside  available  to  the  system 


*I  -  1,ITT;  'J  «•  1,N0TT;  ®K  -  1,NP;  L  -  1,LRT;  -  1, NTOUR;  ^N  »  1,CI0S 


To  input  the  tour  deck  card  B  information  onto  a  permanent  storage 
disk  prior  to  the  simulaL-'.''n3,  the  user  employs  Program  FILL.  Program 
FILL  writes  the  card  B  data  onto  a  disk  in  one  of  two  ways:  it  either 
zeroes  the  entire  matiix  node  and  then  enters  the  non-zero  elements 
into  the  node,  or  it  enters  the  elements  into  an  e:':'*tlng  uatrix  node. 

By  entering  elements  i  'to  an  existing  node,  the  user  holds  ''oiT'pu''''r 
time  at  a  minimum  when  he  is  modifying  small  parts  oi  t..  "ta,  '  i::;  \  6 
illustrates  the  data  deck  setup.  Table  5  describes  the  specif ie  tiauv.  con¬ 
struction  and  Table  4  defines  each  of  the  input  variables  for  Progrr-  TILL. 
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Figure  6.  Data  deck  setup  for  Program  FILl 
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Table  5 

INPUT  TO  program  hl,L 


Card  Description 

No. 

Card 

of 

Var i abl  e 

Columns 

Format 

Dimension 

1 

M 

1-4 

14 

Node  Definition 

1 

ISTART 

1-4 

14 

ISTOP 

5-8 

14 

Nude  Input  Setup  consists 

of  ISTOP- ISTART +1  groups 

,  each  with  one 

Card  A  and  K  Card 

s  B  • 

Card  A 

1 

K 

1-4 

14 

IZERO 

5-8 

14 

Card  B 

K 

N 

1-6 

lo.O 

I 

7-8 

12 

J 

9-10 

12 

PER 

11-12 

F2.2 

1  Repeat  Node  Definition  and 

Node  Input  Setup  cards 

for  each  group  of 

node.s  to  be  input 

to  system. 

Program  Terminator 

n 

-none- 

l-8o 

El  ank 

Table  4 

DEFINITIONS  OF  VARIABLES  FOR  PROGRAM  FILL 

Variable  Definition 

M  Dimension  of  a  square  matrix.  In  the  GMM  data,  this  Is 

equal  to  the  greatest  time  in  the  system.  (maxlen) . 

ISTART  Number  of  starting  matrix 

ISTOP  Number  of  ending  matrix 

K  Number  of  data  cards  to  be  input  to  matrix 

IZERO  If  IZERO  .EQ.  0,  matrix  will  be  zeroed,  then  data  entered. 

If  IZERO  .EQ.  1,  data  will  be  entered  into  existing  matrix. 
N  Value  of  the  matrix  element 

I  Starting  row  of  the  matrix  entered 

J  Starting  columi.  of  the  matrix  entered 

PER  Percent  of  N  entered  into  each  element  of  a  diagonal  vector 

starting  at  element  I,J.  If  PER  .EQ.  0,  lOO^t  enters 
element  I,J. 
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Flows  within  the  system.  Three  types  of  personnel  flow  are  avail¬ 
able  within  the  GMM — direct  movement  between  nodes,  movement  between 
nodes  based  on  node  requirements,  and  movement  to  an  initial  node  position. 

The  initial  and  final  transfer  rules  determine  direct  movement 
between  nodes.  This  movement  can  be  from  any  row  and  column  of  a  losing 
node  to  any  row  and  the  same  column  of  a  gaining  node.  These  data  cards 
specify  which  elements  of  the  losing  and  gaining  nodes  within  the  node 
clusters  and  which  percentages  of  these  elements  will  be  moved  from  and 
to  the  respective  nodes.  Normally,  this  movement  represents  the  Type  1 
flow  of  personnel  across  skill  levels  within  command  elements  prior  to 
the  assignment  process.  For  example,  a  transfer  rule  might  specify  that 
10  percent  of  the  personnel  in  the  eighth  month  of  the  ST  node  must  move 
to  the  next  skill  level  within  the  ST,  or  to  the  outside  of  the  system. 

tiovement  between  nodes  based  on  node  requirements,  described  by  the 
priority-of-f ill  rules,  represents  the  Type  2  flow  or  the  rotation  phe¬ 
nomenon.  These  rules  direct  the  search  and  assignment  procedure  by 
specifying  the  order  in  which  the  program  searches  certain  nodes  for 
personnel  to  fill  other  node  requirements.  For  example,  in  order  to 
fill  the  ST  requirements,  the  rules  might  initially  direct  the  search  to 
personnel  who  have  completed  24  months  in  a  CONUS  node,  after  which 
search  might  be  directed  to  personnel  who  have  completed  24  months  in 
the  LT. 

Movement  to  an  initial  node  position,  i.e.,  to  the  first  row  and 
first  column  of  a  gaining  node,  is  described  by  the  Last  Resort  Tour 
Rules.  This  movement  represents  the  Type  4  flow.  The  Last  Resort  Tour 
Rules  currently  specify  where  unassigned  trainee  input  will  go.  For 
example,  new  trainees  who  are  not  assigned  to  fill  the  ST  requirements 
will  be  assigned  to  a  holding  node  in  CONUS.  Both  data  segment8--status 
of  the  system  and  flows  within  the  system--are  input  once  it  the  start 
of  the  simulation  and  cannot  be  modified.  In  the  SIMPO-IT  GMM,  the  user 
will  be  able  to  modify  these  parameters  during  the  computer  run. 

Special  monthly  system  requirements .  The  last  section  of  input  data 
consists  of  two  cards  for  each  time  period,  one  card  stating  the  node 
quota  or  authorized  strength  and  one  the  input  to  the  system  from  the 
outside.  Based  on  the  inagni:ude  of  these  monthly  node  requirements,  per¬ 
sonnel  flow  through  the  system  according  to  the  priority-of-f ill  rules. 

If  specified  in  the  priority-of-fill  rules,  the  quotas  may  be  partially 
filled  by  the  input  to  the  system  which  represents  new  output  from  the 
training  schools. 

Step  2.  Update  the  system.  System  updating  simulates  the  passage 
of  a  single  time  period  within  the  personnel  system.  Since  both  time  in 
node  and  time  in  system  locate  personnel  within  the  nodes,  the  updating 
procedure  must  include  both  these  time  dimensions.  Therefore,  the  up¬ 
dating  procedure  moves  all  personnel  over  one  row  and  down  one  column  in 
the  node  corresponding  to  completion  in  the  real  system  of  another  month 
in  node  and  in  system. 
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A  reverse  step-wise  process  accomplishes  this  system  updating.  All 
personnel  move  over  one  row  and  down  one  column  starting  with  the  last 
row  and  column  and  working  backwards  to  the  first  row  and  column  in  the 
node.  The  latest  months  In  the  node  and  in  the  system  equal  the  para¬ 
meters  LENGTH  and  MAXLEN,  respectively.  Initially,  all  personnel  who 
have  completed  LENGTH  months  in  the  node  accumulate,  after  loss  rates 
have  been  applied,  in  the  LENGTH  -f  1  row  in  the  node  and  in  a  column 
corresponding  to  their  time  in  the  system  +  1;  personnel  who  have  com-' 
pleted  MAXLEN  months  in  the  system  accumulate  in  the  MAXLEN  +  1  column 
in  the  node  and  in  a  row  corresponding  to  their  time  in  the  node  -h  1. 

The  LENGTH  row  and  MAXLEN  column  in  the  node  are  then  set  equal  to  zero, 
after  which  the  personnel  in  the  LENGTH  -  1  row  and  the  MAXLEN  -  1  column 
are  updated  to -the  LENGTH  row  and  the  MAXLEN  column  of  the  node.  This 
reverse  process  continues  moving  all  personnel  diagonally  within  the  node 
until  the  first  row  and  first  column  equal  zero.  At  this  time  the  up¬ 
dating  process  is  complete. 

Step  3.  Make  initial  transfers .  Initial  personnel  transfers  may 
represent  losses  from  the  system  and/or  movements  between  nodes  within 
the  system.  These  transfers  usually  represent  a  Type  1  flow  across  skill 
or  grade  levels  within  conmand  aggregates.  Since  these  initial  transfers 
occur  outside  the  regular  assignment  process,  they  are  not  a  function  of 
node  requirements.  Instead,  the  amount  of  flow  is  a  direct  function  of 
the  percentages  specified  in  the  transfer  rules.  To  make  these  transfers 
MAINOiM  releases  control  to  subroutine  TIMEUP.  which  transfers  personnel 
from  any  row  and  column  of  a  losing  node  to  any  row  and  the  same 
column  of  a  gaining  node.  Subroutine  TIMEUP  removes  PCT  percent  of  the 
personnel  who  are  in  the  PERU  time  period  from  the  OUTSO  node  within  the 
OUTTO  node  cluster  (see  Table  2  for  definition  of  terms).  These  person¬ 
nel  will  enter  the  PERDTO  time  period  of  the  OUTST  node  within  the  OUTTT 
node  cluster.  In  order  to  specify  movements  representing  losses  out  of 
the  system,  the  gaining  node  coordinates,  OUTTT  and  OUTST,  must  equal 
zero.  Figure  7  illustrates  in  detail  how  TIMEUP  ac:ompli8hes  these 
different  transfers.  In  the  SIMPO  II  GMM,  a  new  subroutine  with  the 
capability  of  transferring  personnel  from  any  row  &id  column  of  a  losing 
node  to  any  other  row  and  column  of  a  gaining  node  will  replace  sub¬ 
routine  TIMEUP. 

To  transform  records  of  the  initial  transfer  movements  into  usable 
form,  Subroutine  TIMEUP  releases  control  to  Entry  Pack  of  Subroutine  X'jQB 
This  routine  PACK,  written  in  COMPASS,  compresses  or  packs  a  cumbersome 
data  description  of  the  transfer  movements  into  a  data  description  which 
can  be  easily  handled  and  stored.  Entry  PACK  packs  eight  fixed-point 
24-bit  computer  words  into  three  fixed-point  24-bit  computer  words 
(Figure  8).  These  three  words  contain  data  describing  the  node  and  node 
cluster  from  and  £o  which  the  personnel  were  moved,  the  number  of  per¬ 
sonnel  input  into  the  node  and  lost  to  the  system,  and  the  time  periods 
in  the  Input  node  and  In  the  system. 
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Figure  7.  Logical  Flow  of  Subroutine  TIMEUP 


WORD  1 


WORD  2 


WORD  3 


Figure  8. 


- 6  BITS  — 

-^6  BITS  — 

- 6  BITS  — 

- 6  BITS  — 

OUTTOUR 

OUTSUB 

INTOUR 

INSUB 

24  BIT  WORD 

17  BITS - -  BUS - 


24  BIT  WORD 

17  BITS  — -  7  bits 


PERSONNEL  LOST  TO  SYSTEM 


TIME 

IN  SYSTEM 


24  BIT  WORD 


Storage  of  data  by  entry  PACK 


-  25  - 


After  PACK  has  compressed  the  data  description  into  three  computer 
words,  TIMEUP  regains  control  and  stores  these  three  words  on  a  temporary 
disk.  This  cycle  of  first  packing  and  then  storing  the  data  occurs  after 
each  transfer  of  personnel. 

Step  4.  Calculate  node  totals.  All  personnel  assets  within  each 
node  and  node  cluster  are  matched  with  node  and  node  cluster  requirements 
in  Step  ^  to  determine  the  number  of  additional  personnel  needed.  To 
calculate  these  asset  totals,  MAINGMM  releases  control  to  subroutine 
SUMMARY. 

Subroutine  SUMMARY  has  the  capability  of  obtaining  sums  of  partial 
and/or  complete  matrices  and  groups  of  matrices.  Elements  In  the  rows 
and  columns  of  any  one  matrix  may  be  added  to  elements  in  the  rows  and 
columns  of  any  other  matrix. 

In  the  GMM,  where  only  sums  of  complete  nodes  are  calculated,  sub¬ 
routine  SUMMARY  parameters  NUMELEifl  and  INDIV  determine  how  these  sums 
will  be  calculated  (Figure  9)<  M/iINGMM  sets  NUMELEM  equal  to  two  and 
INDIV  equal  to  the  node  to  be  summed.  Personnel  In  rows  one  to  LENGTH 
and  In  columns  one  to  MAXLEN  are  summed  for  each  node.  This  sum  repre¬ 
sents  all  personnel  who  are  serving  In  the  node  during  a  particular  time 
period.  It  does  not,  however,  include  personnel  in  row  LENGTH  +  1,  who 
have  completed  LENGTH  time  periods  In  the  node  and  are  waiting  for  new 
assignments. 

Subroutine  SUMMARY  can  also  obtain  sums  of  partial  nodes  or  groups 
of  nodes.  A  MATGRPS  vector,  which  has  NUMELEM  elements,  controls  the 
order  In  which  the  Individual  and  aggregate  node  sums  are  calculated. 

The  value  of  each  element  In  the  MATGRPS  vector  refers  to  the  node  to  be 
suimed.  A  zero  value  for  an  element  directs  the  subroutine  to  obtain  an 
aggregate  sum  of  all  node  totals  calculated  since  the  last  zero  element. 
For  example.  If  NUMELEM  were  equal  to  ^  and  the  values  of  the  three 
MATGRPS  elements  were  1,  2,  and  0,  respectively,  the  subroutine  would 
sum  the  first  and  second  nodes  separately  and  then  add  these  two  sums 
together.  In  this  example,  a  total  of  nodes  1  and  2  could  represent  all 
AUS  and  RA  personnel  serving  in  the  short  combat  tour  area. 

To  obtain  sums  of  partial  nodes,  other  summation  guidelines  of  In¬ 
clusive  rows  and  columns  for  each  node  must  be  Input  as  elements  of  the 
following  vectors:  BEGROW,  ENDROU,  BEGCOL,  and  ENDCOL.  Disregarding 
zero  elements  in  the  MATGRPS  vector,  the  order  of  other  elements  within 
the  four  guidline  vectors  corresponds  to  the  order  of  the  elements  in 
the  MATGRPS  vector.  The  first  element  in  each  of  the  BEGROW,  ENDROW, 
BEGCOL,  and  ENDCOL  vectors  refers  to  the  first  row,  last  row,  first 
column,  and  the  last  column  to  be  summed  In  the  first  MATGRPS  node.  If 
the  first  element  in  BEGROW  were  equal  to  2,  in  ENDROW  to  12,  In  BEGCOL 
to  1,  and  in  ENDCOL  to  48,  then  the  first  node  sum  would  contain  all 
personnel  In  the  second  to  the  twelth  time  periods  of  the  node  and  the 
first  to  the  forty-eighth  time  periods  In  the  system.  This  sum  could 
represent  AUS  non-transients  In  the  short  tour  area. 
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Figure  9.  Flow  chart  for  Subroutine  SUMMARY 


Step  ‘3.  Detennlne  requirements.  To  determine  node  and  node  cluster 
requirements,  MAINGMM  subtracts  the  assets  In  each  node  (calculated  In 
step  4  by  subroutine  SUMMARY)  from  the  node  requirements  Input  for  the 
same  time  period.  If  additional  personnel  are  needed  to  fill  these  re¬ 
quirements,  the  program  continues  to  step  6.  Other,  it  executes  step  7* 

Step  6 .  Fill  requirements.  Once  the  node  cluster  requirements 
have  been  determined,  MAINGMM  relinquishes  control  to  subroutine  FILUP, 
which  searches  for  personnel  to  fill  these  requirements. 

Subroutine  FILUP  uses  an  identical  procedure  to  apply  each  priority- 
of-fill  rule  in  the  order  in  which  each  rule  is  input  to  the  program 
(Figure  lO).  Initially,  FILUP  checks  to  see  if  the  gaining  node  has  any 
requirements.  When  these  requirements  are  positive,  FILUP  determines 
how  the  personnel  search  procedure  will  be  handled.  If  the  losing  node 
is  inside  the  system,  subroutines  FILUP  and  LOCK  search  for  personnel 
within  the  node.  If  the  personnel  are  from  outside  the  system,  sub¬ 
routine  OUTSIDE  handles  the  search. 

When  the  losing  node  Is  within  the  system,  FILUP  first  seeks  per¬ 
sonnel  who  have  completed  LENGTH  tltiie  periods  In  the  losing  node  and  are 
available  for  reassignment.  If  the  requirements  for  the  gaining  node 
are  still  positive  after  using  all  of  these  available  personnel,  FILUP 
begins  to  remove  personnel  before  they  have  completed  LENGTH  months  In 
the  losing  node.  In  order  to  remove  personnel  within  a  node,  FILUP  calls 
subroutine  LOCK,  a  routine  which  locates  personnel  within  a  losing  node, 
removes  them,  and  Inputs  them  Into  a  gaining  node. 

Initially  subroutine  LOCK  searches  for  personnel  within  the  LENGTH 
row  and  column  MAXLEN  to  one  of  the  losing  node.  After  each  group  of 
personnel  has  been  located  and  removed  from  the  losing  node,  subroutine 
LOCK  adds  them  to  the  gaining  node.  In  subroutine  LOCK,  the  time  period 
in  the  system,  or  column,  of  the  gaining  node  is  identical  to  that  of 
the  losing  node.  The  time  period  in  the  node,  or  row,  however,  becomes 
row  one  in  the  gaining  node  regardless  of  what  it  was  in  the  losing  node. 

If  there  are  not  enough  personnel  in  the  LENGTH  time  period  to  fill 
the  node  requirements,  LOCK  searches  in  a  similar  manner  In  the  LENGTH  -  1 
row.  This  reverse  row-by-row  procedure  continues  until  all  node  require¬ 
ments  are  met,  or  until  LOCK  has  searched  in  the  AFTER  row  (Figure  11). 
Control  then  returns  to  FILUP. 

To  decrease  the  amount  of  computer  time  needed  for  searching  the 
entire  node  system  column  dimension  of  MAXLEN  to  one,  LOCK  uses  a 
COMPASS  function  MX. 
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Figure  11  continued 


Function  MX  determines  the  latest  time  period  in  the  system  or  the 
last  column  where  there  are  personnel.  For  each  node  there  is  an  INE 
vector  whose  positions  correspond  to  specific  rows  in  the  node  (MX  sets 
each  INE  element)  with  elements  equal  to  the  last  column  in  that  row 
where  personnel  are  found.  For  example,  the  element  INE(X)  refers  to  the 
X  row  and  the  INE(X)  column  In  row  X  of  the  losing  node.  By  searching 
backwards  in  row  X  from  the  INE(X)  column  instead  of  from  the  MAXLEN 
column,  subroutine  LOCK  reduces  search  time  significantly. 

If  the  losing  node  represents  personnel  input  from  outside  the 
system,  FILUP  relinquishes  control  to  subroutine  OUTSIDE.  Instead  of 
to  subroutine  LOCK.  As  shown  in  Figure  12,  subroutine  OUTSIDE  has  an 
option  of  either  fixed  or  variable  input  from  outside  the  system.  Under 
the  fixed  input  option,  the  amount  of  input  from  the  outside  for  this 
time  period,  specified  in  the  monthly  requirements  Input  card,  is  the 
maximum  which  can  be  allocated  against  all  gaining  node  requirements. 
Under  the  variable  input  option,  however,  the  maximum  amount  of  input 
from  the  outside  is  a  function  of  the  requirements  for  the  individual 
gaining  nodes.  These  two  options  can  either  test  the  operating  level 
of  a  system  with  a  fixed  amount  of  input  or  determine  the  number  of  new 
personnel  needed  to  maintain  a  certain  operating  level. 

These  new  personnel  enter  the  first  row  and  first  coltnim  of  the 
gaining  node,  after  which  control  returns  to  subroutine  FILUP. 

Data  descriptions  of  each  personnel  movement  within  subroutines 
FILUP,  LOCK,  and  OUTSIDE  are  packed,  using  subroutine  X59B-Entry  PACK, 
and  then  are  stored  on  a  temporary  disk.  This  cycle  of  transferring 
control  to  FILUP  at  the  beginning  of  each  priority-of-fill  rule  and  then 
to  LOCK  or  OUTSIDE,  if  necessary,  continues  until  every  priority-of«fill 
rule  has  been  applied. 

Step  7»  Make  final  transfers.  Final  transfers  are  handled  identi¬ 
cally  to  initial  transfers  by  subroutine  TIMEUP.  The  difference  is 
merely  one  of  timing.  The  final  NOFIRST  +  1  to  NOTT  transfers  occur 
after  all  personnel  assignments  have  been  made  against  the  requirements, 
whereas  the  initial  1  to  NOFIRST  transfers  occur  prior  to  the  assignment 
process.  These  final  transfers  round  out  the  system  before  the  next  time 
period  begins  by  reassigning  all  remaining  personnel  who  must  be  moved. 
For  example,  when  a  person  completes  a  combat  tour  (ST),  he  must  be  re¬ 
assigned  regardless  of  the  requirements  in  other  tour  areas.  Thus,  as  a 
final  transfer  movement,  personnel  who  have  completed  the  ST  and  who  have 
not  been  reassigned  might  be  transferred  to  a  CONUS  tour. 

Step  8.  Input  personnel  from  outside  into  the  system.  During  each 
time  period  of  the  simulation,  the  user  may  input  into  the  system  a  fixed 
number  of  new  personnel  representing  new  trainees.  If  this  new  input  is 
not  used  within  the  time  period  to  fill  the  node  requirements,  it  must 
enter  nodes  within  the  system.  These  nodes  are  called  Last  Resort  Tours. 
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Figure  12.  Flow  chart  for  Subroutine  OUTSIDE 


MAINGMM  Inputs  these  unasslgned  new  trainees  into  the  first  row  and  first 
column  of  the  gaining  nodes  specified  in  the  Last  Resort  Tour  Rules. 

This  position  in  the  node  represents  their  first  time  period  in  the  node 
and  their  first  time  period  in  the  system.  As  with  ot-her  movements  of 
personnel,  these  movements  are  packed  and  stored  on  a  temporary  disk. 

Step  9«  Increase  time  period .  At  the  end  of  the  eight  previously 
mentioned  simulation  steps,  the  variable  INTEG,  representing  time  periods, 
is  increased  by  one.  A  new  time  period  within  the  personnel  system  begins. 

After  Step  9^  the  cycle  repeats  itself  for  the  next  time  period. 

The  cycle  ends  when  INTEG  equals  LAST,  the  last  time  period  simulated. 

At  this  time,  MAINGMM  calls  Entry  UNPACK  of  subroutine  X59B,  which 
reverses  the  process  of  PACK  by  unpacking  the  data  descriptions  of  the 
flow  movements  within  the  system.  These  flow  data  and  other  data 
specified  by  the  user  are  then  output  by  an  on-line  printer. 


GMM  SAMPLE  PROBLEM 

A  problem  for  GMM  proposed  by  the  Plans  and  Programs  Office  of  the 
Office  of  Personnel  Operations  concerned  the  feasibility  of  management's 
proposed  policy  limiting  overseas  assignments.  The  proposal  had  been 
made  to  eliminate  all  third  overseas  assignments  except  for  those  persons 
volunteering  for  such  duty  or  for  Individuals  who  had  received  an  inter- 
Jectory  assignment  in  Continental  United  States.  Very  limited  informa¬ 
tion  was  available  on  a  data  base  for  a  sample  subsystem.  In  a  total 
inventory  of  IO78  persons,  3IO  were  assigned  to  short  tour,  II9  to  long 
tour,  and  649  to  CONUS.  Of  the  510  in  short  tour,  100  had  served  in 
long  tour  immediately  prior  to  the  short  tour  assignment  while  210  had 
been  in  CONUS.  The  II9  in  long  tour  were  I07  from  short  tour  and  12 
from  CONUS.  The  system  was  to  be  considered  closed--no  losses  or  gains. 
(Since  an  Intermediate  level  personnel  system  was  being  considered,  this 
assumption  was  not  bad--losses  to  the  real  system  would  be  replaced  with 
persons  having  had  similar  assignment  history  and  the  total  personnel 
Inventory  would  remain  the  same.)  The  question  posed  was:  Will  it  be 
possible  to  maintain  the  authorized  overseas  manning  levels  while  at  the 
same  time  allowing  18  months  In  the  CONUS  tours,  if  no  persons  are  re¬ 
quired  to  serve  three  consecutive  overseas  tours? 

Considering  the  inadequacy  of  the  furnished  data  base,  the  problem 
was  not  complex  enough  for  a  simulation  model.  No  more  Information  was 
available  than  that  required  by  the  steady-state  models  described  in  the 
early  reports  on  manpower  models  by  Sorenson  and  Olson  (5,6).  Pages  44 
and  45  of  the  Nomogram  report  show  estimates  of  the  possible  CONUS  tours 
under  two  conditions:  1)  with, no  sequential  overseas  assignments  and 
2)  with  all  long  tour  assignees  going  to  a  short  tour  assignment  while 
all  short  tour  assignees  return  to  CONUS  for  an  assignment  there.  With 
no  sequential  overseas  assignments,  the  nomogram  shows  that  duration  of 
the  CONUS  tour  would  fall  somewhere  between  I8  and  24  months.  Under  the 
second  condltlon--with  two  sequential  assignments--the  CONUS  tour  would 
be  somewhat  longer  than  24  months,  but  less  than  50  months. 
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Since  SIMPO-I  research  analysts  wished  to  test  the  adaptability  of 
the  G^,  the  decision  was  made  to  set  up  a  model  of  this  problem  with 
the  GMM.  For  modeling  purposes,  the  system  was  structured  into  three 
main  tour  areas,  short  tour,  long  tour,  and  CONUS.  Subtours  were  deter¬ 
mined  by  prior  history;  short  tour  after  CONUS,  short  tour  after  long 
tour,  long  tour  after  CONUS,  long  tour  after  short  tour,  CONUS  after 
short  tour,  CONUS  after  long  tour,  and  CONUS  before  overseas  assignment. 
Since  the  GMM  uses  a  synmetrlcal  arrangement  of  subtours  within  main  tour 
areas,  three  subtcurs  were  needed  in  each  main  tour.  Tour  numbers  were 
assigned  by  the  following  design: 


Previous  Assignment 


ST 

LT 

CONUS 

Present 

ST 

1,1 

1,2 

1,5 

Assignment 

LT 

2,1 

2,2 

2,5 

CONiJS 

5,1 

5,2 

5,5 

Tours  1,1  and  2,2  are 

dummy  tours. 

The  starting  data  show  them  empty 

and  no  priorities  are 

included  for 

transferring  persons  to  them  during 

the  simulation.  Maximum,  minimum. 

and  actual 

starting  tour  durations 

in  months  are  shown  below; 

Maximum 

Minimum 

Actual 

ST 

12 

12 

i2 

LT 

56 

12 

18 

CONUS 

56 

12 

18 

Assignment  priorities 

used  were: 

Into 

From 

Part  of  quota 

Subtour 

Subtour 

After 

to  be  filled 

1,5 

5,5 

56 

100^ 

1,5 

5,2 

56 

II 

1,2 

2,5 

56 

II 

1.5 

5,5 

18 

II 

1,5 

5,2 

18 

II 

1,5 

5,1 

18 

II 

1,5 

5,5 

12 

It 

1,5 

5,2 

12 

It 

1,5 

5,1 

12 

It 

1,2 

2,5 

12 

It 
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Into 

From 

Part  of  quota 

Subtour 

Subtour 

After 

to  be  filled 

2,5 

r 

i  y 

18 

100^ 

2,5 

^,1 

18 

II 

2,1 

1,5 

12 

II 

2,5 

5,5 

12 

II 

2,5 

5,1 

12 

II 

2,5 

5,2 

12 

It 

2,1 

1,5 

12 

M 

2,1 

1,2 

12 

II 

5,1 

1,1 

12 

U 

5,1 

1,5 

12 

It 

5,1 

1,2 

12 

II 

5,2 

2,1 

56 

II 

5,2 

2,5 

56 

II 

For  the  24-month  simulation,  the  manning  levels  for  the  three  ma<n  tour 

areas  were  to  be: 

Short 

Long 

Month 

Tour 

Tour 

CONUS 

1 

510 

119 

67' i 

2 

510 

110 

67' » 

5 

510 

llO 

67'. 

4 

5in 

119 

679 

5 

510 

119 

709 

6 

510 

119 

709 

7 

510 

119 

709 

8 

510 

119 

759 

9 

510 

119 

759 

10 

510 

119 

759 

11 

510 

119 

759 

12-56 

510 

119 

769 

Changes  in  the  designated 

manning  level  for  CONUS  were  made  to  test 

observed  shortages  output 

by  the  model  and  to 

m£.ke  sure  persons  ending 

overseas  tours  would  bf  moved  by  the 

model  to 

a  different  assignment. 

Card  1 

( Parameter  card 
describing  the 
limits  of  the 
sydtem  being 
modeled) 


Data 

Columns 

Deck  Description 

Variable* 

Value  in 

1-5 

( Identification 
of  Sample) 

Sample  Problem 
MAJ 

4-7 

NTOUR 

3 

8-11 

ITT 

9 

12-1‘i 

FIRST 

1 

16-19 

LAST 

20 

20-23 

TFILL 

0 

24-27 

NP 

25 

28-31 

MAXSUB 

5 

32-55 

CIOS 

0 

36-59 

JC 

6 

40-43 

IC 

0 

44-47 

MAXLEN 

48 

52-57 

LRT 

0 

58 

IDISTON 

0 

59 

I PUNCH 

1 

*  See  Table  2  for  definitions 


Cards  2-10  Subtour  Parameter  Cards* 

Card  Number 


Variable 

Cols 

2 

5 

4 

5 

6 

7 

8 

9 

10 

Main  tour 

area 

1-10 

1 

1 

1 

2 

2 

2 

5 

5 

5 

Subtour 

11-20 

1 

2 

5 

1 

2 

5 

1 

2 

5 

Quota 

21-30 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Actual 

31-40 

0 

100 

210 

107 

0 

12 

0 

0 

649 

Nvimber  of 
Subtours 

41-50 

3 

5 

5 

5 

5 

5 

5 

5 

5 

Maximum 
Length  of 
Subtour 

51-60 

12 

12 

12 

56 

56 

36 

47 

^7 

47 

Loss  rate 

61-70 

0 

0 

0 

0 

0 

0 

0 

0 

0 

*  In  this  problem  the 

data  matrices 

were 

already 

on  a 

computer 

disk 

and 

were  called  in  when  needed. 


-  58  - 


Cards  11- 

•35  Priority-of-Fill 

Into  Into 

Rules 

Out  of 

Out  of 

After 

Percent* 

Tour  Area  Subtour 

Tour  Area 

Subtour 

(months) 

Available 

Cols 

1-8 

9-16 

17-24 

25-52 

o 

1 

41-46 

Card  11 

1 

5 

5 

3 

56 

loco 

Nr. 

1 

5 

5 

2 

36 

1000 

15 

1 

2 

2 

5 

56 

1000 

14 

1 

5 

3 

3 

18 

1000 

15 

1 

5 

3 

2 

18 

1000 

16 

1 

3 

3 

1 

18 

1000 

17 

1 

5 

3 

3 

12 

1000 

18 

1 

5 

3 

2 

12 

1000 

19 

1 

5 

3 

1 

12 

1000 

20 

1 

2 

2 

3 

12 

1000 

21 

2 

5 

5 

3 

18 

1000 

22 

2 

5 

3 

1 

18 

1000 

25 

2 

1 

1 

5 

12 

1000 

24 

2 

5 

5 

5 

12 

1000 

25 

2 

5 

5 

1 

12 

1000 

26 

2 

5 

5 

2 

12 

1000 

27 

2 

1 

1 

5 

12 

1000 

28 

2 

1 

1 

2 

12 

1000 

29 

5 

1 

1 

1 

12 

1000 

50 

5 

1 

1 

5 

12 

1000 

51 

5 

1 

1 

2 

12 

1000 

52 

5 

2 

2 

1 

56 

1000 

55 

5 

2 

2 

5 

56 

1000 

*  Read  with  three 

decimal  places. 

\ 
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Cards  34-69  Tour  Area  Quotas 


Tour  Area 


Card  No. 

l(Cols  1-B) 

2  (Cols  9-16) 

(Cols  17 

54 

510 

119 

679 

55 

510 

119 

679 

56 

510 

119 

679 

57 

510 

119 

679 

58 

510 

119 

709 

5q 

510 

119 

709 

40 

510 

119 

709 

41 

310 

319 

709 

42 

510 

119 

709 

45 

510 

119 

709 

44 

510 

119 

709 

4')-0  ) 

510 

119 

769 

Card  70  End  of  file  card. 

To  put  the  data  matrices  on  the  computer  disk,  a  special  program 
called  FILL  was  used.  FILL  develops  a  data  matrl  z  from  Input  consisting 
of  matrix  maximum  dimensions  and  individual  cell  values,  together  with 
the  row  and  column  for  each  cell. 
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Reassignments  actually  simulated  by  the  model  are  shown  below: 


To  Short 

Tour 

To  Long  Tour 

To 

CONUS 

from 

from 

from 

from 

from 

Mont  h 

CONUS 

LT 

CONUS 

ST 

ST 

LT 

1 

50 

50 

? 

50 

50 

■*. 

25 

25 

4 

25 

25 

25 

25 

6 

25 

25 

7 

2“^ 

25 

8 

25 

25 

<) 

25 

25 

in 

25 

25 

11 

25 

25 

12 

25 

25 

I’i 

50 

50 

14 

’0 

50 

15 

h 

25 

16 

25 

25 

17 

25 

25 

18 

25 

25 

19 

25 

5 

25 

5 

20 

25 

6 

25 

6 

21 

25 

6 

25 

6 

22 

25 

6 

25 

6 

25 

25 

6 

25 

6 

24 

25 

6 

25 

6 

2n 

29 

1 

7 

50 

6 

26 

29 

1 

7 

50 

6 

27 

24 

1 

7 

25 

6 

28 

24 

1 

7 

25 

6 

29 

24 

1 

7 

25 

6 

50 

24 

1 

7 

25 

6 

51 

24 

1 

7 

25 

6 

52 

24 

1 

7 

25 

6 

55 

24 

1 

7 

25 

6 

54 

24 

1 

7 

25 

6 

55 

24 

1 

7 

25 

6 

56 

24 

1 

7 

25 

6 
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The  average  CONUS  tour  duration  for  overseas  replacements  was  up  from 
18  months  at  the  start  of  the  simulation  to  over  21  months  at  the  end 
of  the  first  simulated  year  and  to  over  2^  months  at  the  end  of  the 
second  year.  No  sequential  overseas  assignments  were  necessary  during 
that  time.  During  the  third  year  simulated,  the  CONUS  tour  dropped 
somewhat,  but  not  below  18  months  unless  some  persons  were  assigned 
to  two  sequential  overseas  tours.  These  results  are  in  agreement  with 
those  obtained  from  the  Nomograms.  It  is  possible  to  simulate  replace¬ 
ments  to  short  tour  from  long  tour.  Such  runs  were  not  made  with  the 
GMM,  but  results  could  be  expected  to  follow  the  tour  lengths  computed 
for  the  Nomogram  using  the  same  options. 

For  another  application  of  the  GMM  see  the  DISTRO  model  (7)« 
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APPENDIX 


COMPUTER  LISTINGS  OF  6MM  ROUTINES 


Preceding  page  blank 
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^  ^  ri  n 


3I/l?/130o  fortran  (3  D/mSOS 


PROPpaW  MAINGMM 

APRIL  iRft?  N'ARVA 


FII.UP 

NARVA 

April  1R67 

OMTSlOF 

NARVA 

April  1907 

LOCK 

NARVA 

April  1‘’0  7 

COMMON 

NOFTPST  «NOl.AST 

♦  Not  T 

.ITRtIME  ♦integ 

COMMON 

IP  .INP(48) 

»  IEMD 

.OUTtT (200) 

COMMON 

OllTST  (POO)  ,OuTTo(200) 

.OUTsO(20n) 

COMMON 

PCT(?00)  .PERDIPOO) 

♦  J»^ 

,SYST(4ft'.A8) 

COMMON 

PA(IOO)  ,NaV(100) 

♦Length) 1 00)  ♦out(Ioo) 

COMMON 

lMTnUR(??0) 

♦1nSUR(220)  ♦OUTTOuR(220) 

COMMON 

lOfMPPO)  .OuTSUH(2?0) 

.AFTcR(220) 

COMMON 

PFH(?P0)  ,RaTE(I00) 

♦  NEEoSdo)  ♦NE(IO.IO) 

COMMON 

lOSOO)  .NEDOO) 

♦NeE(IOO) 

♦IFIlL  ♦MAXLEn 

COMMON 

NTOUR  ,NP 

•  ClOS 

.PERdTO(20o) 

COMMON 

Pnw  «LSTRSTn(10) 

♦LSTnSTTf 10) 

COMMON 

LSTRSTSdO) 

♦  ITT 

♦  HAXsUB  ♦ACTdo) 

COMMON 

NSIT(lO)  ,INOS(10»10) 

. IGRaDEIjZo) 

COMMON 

REP(??0)  .1TYPE(220) 

.PERoHT(220) 

COMMON 

lOISTON  ♦ISUM 

♦GRPSUM(ion)  *PR10(iOO> 

COMMON 

HFOnUi^  ( j  no) 

♦EndROW(IOo)  ♦BEGCOl(IoO) 

COMMON 

FNOCOL ( 100) 

♦  M/^TSUM(10 

o)  ♦MATGRpSI 100) 

COMMON 

Tyre  '  1  oo)  ««:um  (100) 

♦NuM(lOO) 

♦  ACTjjALdOo) 

COMMON 

NPRlKV  «NT 

♦  Ihoi  d 

♦LEN  ♦LEVEL  ♦M 

COMMON 

NCRNOOF ( 1 00) 

♦MtN( 100) 

.GRPtNPRdoO) 

COMMON 

MAXDEPL(IOO) 

COMMON  THt ( 1 000) 

TNTFCiFR 

TYPE  ♦SUH 

♦Actual 

♦  ENOcOl. 

InTFGFR 

CTOS  ,SyST 

♦OuTSUH 

.outtour  ♦after 

TnTFGFR 

ACT  ,0UTT0 

♦OuTSO 

♦outtt  ♦outst 

Integer 

PEOO  ,G0PSUM 

♦BeGROW 

♦ENDrOW  ♦BEOCOl 

InTFGFR 

FIRST  ,PEROTO 

definitions  of  tfpms  useO  in 

1  StmULATTOM  program 

r 

r 

r 

r  T0=  RUN  TnFNTlFlCATlON 

r 

C  NTOU«=  number  of  overall  TOUR  TyPES 

r 

C  ITT*  NIJM0F.R  OF  TOURSITOTAL  NUMOj-p)  TO  RE  SIMULATED 

C. 

r  FIRSTrBFGiNNTNG  TIME  PERIOD 

r 

c  haxsub=maxtmiim  numrer  of  subtours 

r 

r  NPs  total  nmmhfr  of  priorities  input 

c 

r  NOUTINTOIIR  X  MAXSUB)  =  total  output  from  each  tour  and  StWTOUR 

C  AVATLABLE  FOR  ANOTHER  TOUR 

r 

r  LOSSES  (NTOUR  X  MAXStlB)  *  TOTAL  NUMBER  OuTPUt  FROM  £ACH  fOUR, 

c  subtour  lost  from  system 

r 

C  HIDOUTI  NTOUR  X  MAXSUB)  ■  TOTAL  OUTPUT  PROM  EACH  T0UR»  SUBTOUR 

c  OURTNG  month  available  for  another  tour 

c 

t.  MTDLOSSI  INTOUH  X  MAXSUB)  r  TOT^L  OUTPUT  FRqM  FaCH  TOUR,  SUBTOur 

c  ourtng  month  lost  from  system 

r 

Preceding  page  Wank 


r-  If  II  I  =  PARAMFTEM  tiSFn  TO  OFTFRmInF  IF  qUOTaS  *rE  FOP  TOUR  OR 

r  TF=  n  H»vF  NTOIIW  OUOTAS 

r  TF  lFSS  THFN  or  ORFATFR  THAN  a  HAVF  tTT  qUOTaS 

r 

r 

r  ri«wiiMRFR  OF  CRITTFAI  tours 

r 

r  r?s  NMMHFW  OF  CRITICAI  TOUR  TyPeS  *  ? 

r 

r  nos®  NUMHFR  OF  catfoorif.s  T^JPUT  from  outsIdE 

r  MOWS 

r 

r  maxlfnsmaxtmum  length  of  timf  im  system 

r 

C  Is  J  ,TTT 

r 

r  TYPF(i)rTOiiR  Type 

r 

r  sim( T)  =si)»<ro'iR  typf 

r 

r  OUOTA(l)sOUOTA  OF  TOUw  OR  SUHTOur (DFPFNdS  On  IFILL  PARAMETER) 

r 

r  ACT'IAI  (I)=ArTl)AL  MIMHFR  OF  MFN  in  tour  oP  SuBTOuRF 

r  OFPFNf'FMT  UPON  IFII  \  parameter 

r 

r  NUM(  n  =NUM9i-:j  OF  SUBTO'IHS  !N  TOur 

r 

r  LFNOTH(f)=  LFNGTH  of  S'IHTOIIR 

r 

r  OUT(f )=PERCFMTAGF  output  form  System  after  \  YEaP  EROM  SURTOUr 

r 

r  PeROuTs  PtwrFNTAOF  OF  nuTOUR  AVAIL  FOR  FILL  TO  INTOUP  NEEDS 

r. 

r.— —  ttype  s  type  oe  movement  into  new  Tour 

C— — 0.1  •?  »  (1.1)  INPUT 
r— —  3  «  (1,1)  Input 

r, - 4  r  (1,1)  Input 

r 

r 

Ml cNPsNlsO 

READ  1 ,  in.NTOUR.rTT.FTWST.LAST.IFlLL.NP.MAXSLiR.JC.IC.CIOS.MAXLEN 


)  .LP^  .POW.iniSTON,  IPIJNCH 

C  TF  IFTLL®) .Quotas  are  for  NOOES.IF  IFILLsO.uDOTaS  ARF  for  node  CLUStFrS. 
IF  (IFTLI  -1 )53,54 
S3  NT*NT0UR 

GO  10  55 
S*  NTsITT 
55  roNTTNDE 

c  peaos  In  distribution  parameters 

lF(InTSTON)2l4n,2l40,?l4q 
?14R  MaTOUT«MAT IN=1 

call  input (MaTOUT, MATIN) 

2140  no  2104  IsI.maxlen 

P1A4  IMF(I)»0 

no  2105  1  =  1. ITT 

2105  call  RANWRITF( 11 , INE.MAXLEN.I) 
month*  first 

.)K  =  MAXI.FM*MAXLFN 
00  3210  1  =  1. ITT 
r 

c—— matrix  data  rfao  in. 

REAO  4,  TYPE(T)  .SUR(I)  .lOlJOTE  ,  ACryAL  ( I )  ,NUM  (  t  )  .LENGTH  ( I )  .OUT  ( I ) 
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IJK»  TVPF ( I ) 

MSITd  IK)=NIIM(T) 

C - Tr  OM  (EQUAL  TO  D  RUN  wll.L  RESTART— USE  COftPUTgR  PuNCHEn  MATRIX  0/itA, 

IF  (SswTrnr (4) .EO. 1 >  Ro  TO  1210 
GO  TO  (HOOtROl ) .SSWTCHF O) 

BOO  no  ?  K=i .maxlfn 

2  read  3.  (SYST(  I.K) «Js1 .MAXLFN) 

RQ  TO  802 

801  CALI  RANBEADd^.SYST.JK.I) 

802  CAIL  RANWRdEdO.SYST.  IK.d 
32)0  rOMTlNIlF 

- TpaNSFEP  PRIORITTFR  RFAn  IN. 

read  3. MOTT.MOFIP'^I  .NOLAST 
PpTNT  70.MOTT 'Nrir  I  WST  .NOLAST 

mottsniimbfr  of  possibif  transfers  bftwfen  mOS 

mOF TRSTsNIIMrER  of  transfers  PRIOR  TO  OTHFR  ASSIGNMENTS, 

mOLASTsMUMHER  RF  transfers  following  other  assignments, 

TF  OHTTI  OR  TOMD  TO  WHICH  PERSONNEL  ARF  SEnT  =  o'.  RFmOVALS=LOSSFS  TO  SyStEM 
00  ?n 3  1  =  1  .nott 

REAP  4*  OUdTd  )  .OUTST  d  )  .OlITTOd  )  ,OUTSO(  I)  .PERod)  >ERUTO  ( I )  ♦PCT 
Id) 

2113  PRINT  71.  PCTI  I)  .OlITTOd  )  .OUTSOd  )  .pEROd)  .oUTTt  (1)  .OUTSTd)  , 

IPEROTOd) 
prTnT  7,  1(1 

- priorities  read  in. 

00  5  Isl.NP 

head  6,  INIOUR ( I ) , IMSUH ( 1 ) .OMTTOUR , ? ) .OUTSuR ( !> . AFTER ( 1 ) ,PFR(T ) , 

1  PFPOUT(T),  RFPd).  lOOd).  IGHAnE(I).  ITYPtll) 

NPFPsPFRd)*10n 

print  R,  NPFR  .OUTTOlIRd  )  .OtJTSUBd)  .AFTFRd)  «  INTODRd  )  .  INSIIBII' » 
lITYPFd) 

5  CONTINUE 

- LAST  RESnnT  TOURS  READ  TN, 

on  500  1  =  1  «l  RT 

READ  515.  (LSlRSTOdl.LSTRSTIdl.LSTRSTSdM 
500  PRINT  20.  (LSTRSTO(  I)  ,LSTR‘;TTd)  .LSTrSTSII)  ) 

percents  vector  of  percentages  to  be  APpLIEo  TO  CALCULATED  INPjT 

- PPTNT  SVSTFM  .AS  IS. 

print  1006 

no  807  T  =  1 . 1 T  r 

CALL  RANREADdo  .SYST.JK.T) 

Cali  fltp 
print  135S.  T 
TaT  =  0 
JT  1  =  0 

no  800  Jl  =  l IK, maxlfn 
JTl  =  Jd*1 
I  T=JI»max|  EN-1 
no  808  JrJI.IT 
If(5YST(,U)  009. 808. BOR 

809  IAI=IAI.3 
lBTdAT-2)=  )-JT*l 
TBTdAI)=5YST(J) 

IBT  (lAT-l )=  JII 

808  CONTINUE 

TFdAT)  010. 810.811 
811  print  604.  (101  (,l)  ,  J=1  .TAd 

810  CONTTMUE 


4Q 


r>  ^  D  o 


807  rOMlTMlip 


r-----0rf;lN  MONTH  0Y  MONTH  stmulation, 

00  401  TNTFr,=MONTH«LAST 
TP=n 

DO  372  T  =  1.NT0t)0 
ftCT ( I ) =0 
MEEU<;(T>s<’ 

372  rONTlNUF 

r 

r  UPIiATF  FnTThE  system 

pptnt  ip.tntfo 

00  in  Isl  ,TTT 
Mr  ( T ) =n 

Cali  ranrfao(io.syst,jk,h 
0  =  1  FnOTH(T) 

PPsOlIT  (  T  )  »0/l  ?,0 
LAst  fnotm<  1  )  «] 

TF  (I  A.OT.MAXLEN)  I  A  =  MAX|>:n 

I  FNNOnF  =  |  F^I(iTH  (  T  ) 

TF  (I  FMNOliF.fiT.MAXLFN)  LFNN00E=MA XLEn 

r 

r  UPOATFS  PFRsnNN«^L  MATRICES  BY  COLLECTING  AlL  WmO  HaVE  COMPLETED  ThE  nOqE 

r  In  THF,  la  MOoF  TIMF  MFRIOd  AMO  The  1L»1  system  TTME  period. 

r  RFmOVE  rOMPiETEO  NOOF  PERSONNEL  To  SYSI (L A*  IL ♦  1 ) 

00  14  IL=T«MAXLEN 

II  suM=n 

00  SB  K=l  FNMnOE  .'•IAXLEN 
IF  (SYST  (K,  TL)  )  SF,,SH,S7 
SF  SvST(K,ILI=« 

GO  lO  50 

S7  X1=SYST(K,IL)»PR 
YrSYST (K,TL)-X1 
T.lsXl 

IF(X1-1, 1)11,11,1? 

T?  iKsY 
YsTt<*T 
X1=I.I 

11  iFd.FO.S)  print  4?5,T.t<,IL»SYST<K,Tl:) 

SYST(K,TL)=n 
Tl  SUMS iLSUMtY 
rONTlN'lF 

iFdL.GT.l)  SYST  (LA,  Tl)  sIPFTAIN 
IrFTAINsILSIJM 

all  OTHEPS  ARF  moved  IIP  ONE  MONTH  TN  THF  NODE  aND  ONt  MONTH  IN  THE  SY5>.eM 
kolo  and  iloidspow  and  column  PRioq  TO  updaTIn^  System, 
knew  AMD  TLMtWsROW  AND  COLUMN  AFTEr  UPOATInG  SySTFM. 

ILOLOsMAXLF.N-Il 
TLMFmsMAXI  FN-1l*1 
no  14  ks),i.ennoof 

KNFWsI  FNN0DE-K,1 
KOLOsLFMNOOE-K 
TF(KNFW,I F,l)SR,Sn 
SR  SyST (KNEW. TLNFN) =0 
GO  TO  14 

so  IF  dLNFW.LE.l  IGl  ,G^ 
si  SyST (KNEW, ILNtW) sO 
GO  TO  14 

F?  SYST (KNEW.ILNFW) sSYST (KOLO. ILOLO) ♦SyST (KNEW, ILNeW) 

SyST (KOLD. ILOI  0) =n 
)4  cD'"TTWUE 
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no  n  n  n  n  n  o  n  r»  r»  n 


SY«;i  ( 1 ,1)  =n 

CALL  HAMk<wlTr(in.SYST*JK«T» 

10  cnNTiNiir 

TIMFIIP  accomplishes  TQANSEEhS.  THfSE  transfers  Are  effective  Prior 
Tn  APPI  tcatton  of  fill  oulfs 

TTPTTMFsTRANSFFP  TImF.  IFsi*  all  MOFIRST  transfers  are  made. 

IF=?«  MOFTRST*!  to  MOTT  TRANSFERS  ARF  MAoE. 

ITPTTMFsI 
CALL  TTMEHP 

SUBROUTINE  SUMMARY  CALCULATES  ACTUAL  FOP  EaCH  j^ATRIX 
no  f'k  1=1  .  ITT 
1 ACTUALsACTUAL (  I  ) =0 
call  SIIMMARY(?,I,1ACTUAL) 

ACTUAI  (T)=TACTUAL 
KKKeTYPFU  ) 

ArT (KKK) = act (KKK) ♦actual ( I ) 


yPOATFS  IMF  vector 
CALL  RANRFAn(  11 INE.MAXLEN*  1) 
no  El  K=1 .MAXLEN 
TwE(K»*0 

no  A1  J=1  ,MAY(  FN 

If  (SYST(K,.i)  .GT.o)  INE  (•<)=J 

63  COMiTMUF 

r Al  L  PAMWRITE (11. INF.mAXLEN. I I 

64  roMllNlIF 


26 


27 


PPTMT  ?4 

MAXMOOFrMTMMO()E  =  o 
no  ?7  IsI.NTOUR 
no  ?6  itsi.maxsuo 
MAXMOnF  =  MAXMOIJt*  1 
MINMOnFsMAXSllPo  (  I-l )  ♦  1 
PRINT  ?O.I.ACT<I)  .  (U.  actual  (ID 


A  llOA> 


IlrMlNNORF.MAXNOOE' 


CAlCUl  ATE  minimum  personnel  WITmTN  tour  AREA. 
TF ( IOTSTOM) ?1 4?,  pi  42, pi  43 
21A3  MaTOUT=MAT IN=1 

CAll  ISUMAW(maTOUT, MATIN) 

C 

C—— Input  pfrsonnfi  pfouirements. 

21A2  lF(lFTLI-n  Il.TP.lP 

3)  REAP  33.  (NFFOS  (I  )  ,L=  1  .NTOUR) 

Print  4fl.  (nffus (I  )  ,l  =  i  .NT) 

r,0  TO  14 

32  PFAU  13,  (MF  (I  )  ,1  =1 , 1  tT) 

TS=n 
mO=  1 

no  l^  L=I .MTOUR 
mX  =  NSIT(I  ) 

TS=TS.mx 

no  1*,  K=Mn,TS 

36  MFFOS  (L  )  =''KFnS  (L  )  ♦NE  (K  ) 
wO=TS+ 1 
3b  COMTIniif 
34  CONTTNUF 

r 

r - Tmput  nuT<^li)E  oata 

TF(CTns-n  610.^11.611 
61)  Rlau  11.  ( 1061 n . T  =  ) .rins) 


1  - 


cue  c  u 


f>l  0  TofgT  InuF" 

no  1  fiS  1  =  1  «mToii(j 

ME:ri(T)=MtEDS(T) 

mFF  0«;  (  T  )  =NFpn«?  (  T  1 -AC  T  (  T  ) 

TFlNFFOSd))  17^.16S.lf.S 
MPFOSlIlsn 
IAS  rOMTIMMF 

pnIMT  4q,  (NFn(  I )  ,  1  =  1  .NTO'IQ) 

PdTMT  10,  (NFFf)S(T)  »I  =  1  .NrniiR) 

□0  ISA  I=1.TTT 
mEE ( I »  =Nr ( 1 ) 

►IP  (DsNFdl-ACT'IALd) 

TF(N^d))  177,lAh»lAA 
177  MPd)=n 
ISA  COMTIMIIF 

PpTMT  40. <NEE d  1  . 1  =  1 , T TT) 
pr>lNT  41  ,  INE  d  )  ,  1  =  1  ,  ITT) 

r 

r  COMlpoiS  MFFnS  TO  rtf  Fld.F')  ACCOHOImO  TO  PRIORI  tY  Op  Fill.  LEVEL*! 

PRIM  s?.ini«;Too 
IP  dOISTOM)  ?144,PJ144,?1  4S 
PUS  no  ?U4  l.F vFl  =1  .MPRI.EV 
MATOllT  =  MflTIN=l 

Call  mooii y (matoiit»matim) 

PI  44  CAI  I  FIl  IIP 

TPAnSFPW  all  UMASSIGHEO  personnel 
I TRT IME=P 
CAI  L  TIHEUR 

TEirins)  aip.au.ah 
A1 3  no  SlO  1=1 ,LR1 
Isl  STRSTOd  ) 

K=(  STRSTTd) 

I  I  =LSTPSTS( 1) 

TEdOSlJl)  SlO.SlOtSPO 
S20  Ir=  (K-1 1  *MAXSIJR*1  L 

TALL  RANHEAO  (1 O.SyST , JK, ICl 
SYSTd.nrSVsTd  ,  n  ♦ins(j) 

IP=IP*1 

PRINT  loos,  K,LL , J* ins ( j) , IP 
CALL  PaCKIK.LI  ,0, J, inS(J) ,0*1,1 ,N1.NP,N1) 

CALL  RANwRiTFl  n,N’.  ,1,  IP) 

CALL  RANwRlTFdO,SYST,JK,lC) 
actijai  dr)=ACTiiALdC)  ♦ins(J) 

ArT(K)=ACT<K)*10S(J) 

insi j) =0 

SlO  roNTINUP 
61?  CONTINUE 

- Calculate  oeployarle  personnel  for  pach  tour  arpa, 

ir (IOISTON.LE.O)  go  TO  17 
MATOijT  =  MATIN=l 
CALL  TS'IMAR(MAT0IIT*MATIM) 

17  CONTINUE 
C 

PRINT  1S6 

print  1006 
no  3S0  1=1, ITT 
print  liss*  1 

tall  RANREAni 1 O.SYST , JK , I ) 

Call  flip 


n  o 


I  A  I  =0 
JT  T=o 

no  C,nO  11  =  1  .  IK.MAXLEN 
IT  I =  II  T  *  1 
n=  IT ♦•'lAxi  en-1 
no  #,nn  i=Jl,n 
IK  (  «;ysT  ( .1)  )  (Sni  .f,00,fS01 
KOI  TaI=TAl*l 

lai  (TAI-?)=  .I-J1*1 
I0T  (TAT  )  r'lYST  (J) 
laT  (TAT-1)=  JIT 
KOO  roMlTMIIE 

TKITATl  6O2,K02.K0l 
K'l  PmIMt  K04,  (  TmI  (  I)  ,  J=1  .  T  at  ) 

Ko?  rnnTimiE 

I'-w)  roMTTMiie 

no  T  =  1  .ITT 

TAIL  WAMI^T  AT)  (  1  n.SYSI  ,  JK  4  I  ) 
no  H04  K  =  I  .MAXI  EO 
(MF (K) -n 

no  R04  .lrl»MAXL£N 
TF  (*;y«;t  (K,, in  b04*HA44  8i'i5 
00b  I^'F  (K)  =  I 

B  C  A  roMT  I  'll  ir 

faOb  ("all  PAMWI'I  IF  ( 1 1  .  TNE.MAXLFM,  1  ) 
PPIMT  on-? 


DO  «)no  .TH 

CAT  L  PAMDT  An(  1 1,M1  »  T..I) 

CAT  L  OMF'ArK  (Ml  ,N?,n1) 

Qon  PPTNT  001.  (Ml 


POINT  TSK 

GO  TO  (1?1  1.401)  SSWTrHF(51 
3211  woTTF  (bP.T?l2) 

PAIISK 

IF  (SBwrri-iMb)  .10,1)  r.n  to  4o? 

401  rn'iTT'^illF 
4  02  POT*  T  4  4  4.  i'll  KG 

Tt^  (  1  Pii'jrn.l  I-  .0)  GO  TO  KG 
DO  b  1  1  -  1  .  I  TT 

bl  Pij'if  ■(  *•■  .  (VOT(i).  SliU(T).  IODOTF 
1  O'll  (  (  ) 


.  ACTIJAL(I)»  NIIM'T).  LENgTH(T). 


r 

r  f  fll  1  G  ■  imM' ml  (f|F  AI.I.OCATF  WHTCH  actually  OTSTPtRUTES  PFR«inMNK|,  wTthT 

c  Toiip  AOPAS  lO  COMMAMn  flements 

KS  TF  (  TnTSV'>f')214G.>>14K.?l47 
?147  MA  TOl  I  T=MA  I  TM=  1 

Tali  T  Ann  (MATOIIT. MATIN') 


Xt,OA  KOPMATS 

I  Kno'-' A  I  (  A  1 .  1  2  T  4  ,KG.4,?I  1  1 
1  KnOMAT  ( ?n  1  4) 

4  FORM  A  r  (  K  T  1  0  .Kf,  .  -1) 

6  Format  (MH.  'ifk.I.ih  > 

f  FORMAT  ( TPWIF  T|  I  PRTORTTTKS  FOR  S T MUl  A T I ON , Aa // ) 

8  format  (  TX  .  I  T.  1  x.7HPEOrFNT,SX»F>H  FROm  .T?»1H..I?.SX,tH  AFTER  Tl,7HP 
UOTO'IS.GX.GM  I'lTO  ,  IP.  IH.  ♦  I2»*iX,  15H  MOVFMFNT  TYpE  TP' 

18  forma  I  ( //I  X47H  fMilMMFR  AVATLARLF  for  aNOTHFR  TODr  AFTPP  MONTH.  I4) 

,’0  FORMA  r  ( 1  x,4bhrf  .iain()fr  of  systfm  input  from  outsIof  category  it.sh 

1  INT  0,  (?.  I M. .  I ?) 

?  l  format  . //;»«iiohfosonnFL  TNOICATOR  VECTORS  /Sv » 34i-iVAnjES  =  MOmTH  IN 
ISySTFm  or  r('MJ'T'l/SX.l?HPOSlTIONS  =  MONTH  TN  TOUr  OR  ’ROW) 

24  Format  ( //29'(OflOUf  CLUSTER  ANo  NOOF  TOTALS) 
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2'i  Format (/sx.nHNnHF  Cluster  tiStiH  *  , i6/(13x»5hnode  tizoH  s  »i6)) 
30  format  (//?SH0NnnF  cluster  shortages  x.’iOIlO) 

33  FORMATdOTB) 

40  format ( //I AHONODF  QUOTAS  = » 3K . ( 10  1 1 o d 

41  format ( //I THONOOF  SHORTAGES  b*(10I1o)) 

4B  format  ( //??H0N0DF  CLUSTER  QUOTAS  =»3X.l0n0l 
50  FORMAT  (6I10»  F6.3) 

52  format ( IlHOlnlSTON  =  ,11) 

70  format < IOHOTOTaL  Tf  \NSFEHS  *  ,11,2X.?1H  INITIAL  TRANSFERS  =  »Tl,?X 
I.IRHFINAL  transfers  =  ,11) 

71  FORMAT(1X,Fft.3,(SH  FROM  1 1 ,  I H  ,  ,  1 1  ,2X  .  feHAFTER  *I3'.I4H  TIME  PERIODS  , 
IIhTO  ,n  ,ih,,Ii,2X,AHAFTER  I1,13H  TtmF  PERIODS) 

356  FORMAT!///) 

404  format  (23H  LAST  PERIOD  STMULATE0=  l4) 

425  FORMATlflHoMATHIX  ,T2,3X,4HR0W  , 12, 3» , THCOLUhN  ,t2*3X,2H  r,l4) 

515  F  »RMaT(1I5) 

604  format  (7(4X,?U, 17)  ) 

901  FORMAT(«I10) 

902  format  (//4X6HT0URlN,5X5HSURlN,3X7HToiiR0UT,4X5HSuR0UT,2XnHN  mEN  IN, 
14X6HN  LnST,lX9HM  IN  SYST,1X16H  M  IN  TOUR  (ROW)  ) 

1005  format  (1?H  main  2lfl,BX,4T8) 

1006  format <//l X,  15HPERS0NNEL  NOOES/H022h  SYS  TOyR (COL »ROW) ) 

1355  F0RMaT(1X.I4) 

3212  format  (103H  FOR  ANOTHER  PASS— TURn  OFF  SWITCH  5  AmD  PRESS  MI  / 

1  mT.  to  terminate  press  MT  /  MI, 

2146  STOP 
END 

FORTRAN  DIAGNOSTIC  RESULTS  FOR  MaINGhM 


NULL  STATEMENT  NUMBERS 
23 
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H  o 


Sl'RRniJTTNF 
-Ai  L  rnwMDN 
COMMON 
COMMON 
COMMON 

common 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

common 

common 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

common 

TmTFGPR 
TnTF  gfr 

INTFGFR 
TmTFGFR 
IMTFGFR 


■^l/3?/130r«  fortran  (3.1)/mSOS 
LOCK  (I»J.»'nN.TLOsS*!OAlN.IOUT.lFF,Mix«IrC«lE,Kl,K3  ) 

'’^NOCTR^T*’  IDENTICAL  WITH  ftTHER  SUBROUTINES 

NOFTRST  , NOLAST  »NoTT  .ITRtJMF  flNTEG 

.lNF(4fl)  «lrMn  .OUTtT(?oO). 

oiitst  (pnn>  ,oiiTTn(,?0f))  .0UTs0(?nn) 

RCT(?.)0)  .RFRDlPno.Jx  .SYST(4B.4fl) 


Pa(lOO)  *NA\/(iaO) 
lNTnijR(??n) 


»LrMGTH( 1  On) 
«ImSUR(??0) 


Ion(2?o)  ,0|(TSUB(2?6) 
PFP(220)  .RaTEIIOO) 
lOS(lO)  «NEO(iO)  «NfF(100) 
NTODR  ,Nr‘  .CTOS 

Pow  .LSTRSTO(IO) 

LSTRSTS(IO)  ,I7T 

USTTdO)  .INOS(IO.IO) 

RFP(2?0)  .TTYPF(?20) 
lOlSTON  ,ISUM 
BEGPOW(inO) 

enocolcioo) 

Type  ( loo) .SDR  <100) 

NPRi.EV  ,NT 

NCRNOOE (100) 

MAXOEPL(IOO) 


.OUT  ( 1  (iO) 
.OUTTOuR(??0) 


.AFTrRlp^O) 

.NEEpS(ln)  .NEno',10) 
•IEIlL  .maxlen 

.PERpTO(700) 

.LSTrSTT'(10) 

.MAX^UR  .ACT  do) 

,I6RaOE(?20) 

.PERoUT(;>?0) 
.6opSUM(100)  .PRiOdOO) 

.EmoROwcIOo)  .REGCOl ( 100) 

.MATSUMdOn)  .MATGRPS(ioO) 

.NuMdOO)  .ACT^ALClOo) 

,len  .level 

.GRPiNPRrloO) 


.  IhOLO 
.MiNdOO) 


•  M 


14 


r 

r-— . 


4! 

4  0 
67 


IS? 


r-7-  . 

; 


type 

CIOS 
ACT 
PERU 
PEROTO 
OKK 

IF  ( TGRaoe ( IgNo)  .LE.O  ) 

Print  ions,  iGRAOEdENO). 

Call  trnsposf 

CONTINDF 

Nl=N2sN3=0 

■MON  r  TIME  period  after  WHICH 

mOMTHsMOM. 1 

tttso 

L* ( T-1 ) WMAX. ( 


,SDH 

,syst 

.OUTTO 

.GRPSUM 


GO  TO 
IEND 


.Actual 

.OuTSUB 

.OuTSO 

.BeGROW 


.ENDcOL 

.OUTtOUR 

.OUTtT 

.ENDrOW 


.AFTER 

.OUTST 

.BEGCOl 


14 


removal  takes  place  (maxlen  to  MON.Ii 


.10(0=  TF  ON,  REMOVAL  IS  FROM  PERIOD  FOLLOWING  MOn  ONLY 
iFdOOdE)  -11  40,41,41 
•Ms  ROW  OF  REMOVAL  FROR  THE  MATRIX 
MsMONTH 
GO  TO  47 
MsLENGTHd  ) 

continue  -  -•  --  . 

LENsMAXLEN 

IF(m-monthi  I. ISP, is? 

Tsu=0 

41  !  ‘^4NRe4n  M  >.  INE, MAXLEN. L) 


r.  c 


ISR 


'F  COl.  s  WITH  ♦  Value  In  aLL  rows  GREATFR  than  Or  FO'JaL 

'  r  F  , 

' ’ 1 S4 . 1 S4 

■It'. 

•  .  VAtMF*  *  (TOUR)  .  AluE  OF  iNE(x)>  HIGHEST  COL  oE  RoW 

-Kii-  dV’  entity  tn  or  beyond  removal  rowT - 

Is  dNE(cl))  158.150 
If  ((LI)-mONTH)  150. 151. 151 
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isi  i'sU=i‘;iio 
rONTTMlIF 

C1-  — -FUNCTION  MX  nFTFRMiNFS  hIOHEST  COL  TM  MATRJv  WHtCH  hAS  A  POSITIVE  v/aLuE 
Cl— —IF  loo  1*;  OFF. 

MMN=MX ( INE.MAXLEN) 
r,n  TO  is^ 

CO - ChFCK  THF  NliMhFP  OF  COL  S  WITH  ♦  VA(,UE  IN  AlL  RqWS  EOUAL  TO  MONTH  if  iOoON 

1^4  no  ISFi  11  =  1  .maxlfn 
If  (iNFlLin  1SQ.1*S6 

l^R  IF  ( (I.  I)-MONTM)  iSetl'^T.lSO 
1=>7  1SU»TSII*1 

leif.  continue 

C - Nl^Mr  starting  POW  of  RFmOVAL  FROM  TmF  MATRIX  (Nn/tN  TO  COLl  OF  ROW  M-rFmVlPT 

NmN  =  INF (M) 

15b  lF(lSU-n  1.10?. 10? 

c 

10?  LFNrNMN 
NMN>0 
MMl FNrO 
1PR=1 

C4— SfARCR  for  highest  position  OF  A  POSITIVE  InE  VaLUE 

C4 - NmLEN  is  TMF  HIGHEST  POSITION 

C 

00  ^01  »LKrMONTH,MAXLEN 
iFdNFdl.K)  ,LF.  0)  GO  TO  101 
1F(LFN-TNF  (lLK))101.in?«ln4 
in?  Nmn=TLK 
GO  TO  301 

104  NmLFN=ILK 
301  CONTINUE 

1PRS|PR. 1 

tfinmn.fo.o  .ano,  nmlen.fo.o)  go  to  1 

TF(NMN)  103.301.304 

303  lF(lOn-l ) 104,1,1 

304  1F(NMN-NMLEN)307.305.10S 
107  NmN*NMLFN 

105  MsNMN 

GO  TO  3?? 

C 

121  CALL  PANRFAOdO.SYST,  JK.L) 

CALL  RANRFAO d 1 , INE, MAXLFN. L) 

1F(1nF(M)-1 )q.l?4.1?4 
324  lFN=TNE(M) 

TSUsLEN 
GO  TO  l?l 
C 

12?  CONTINUE 

-  ?  CALL  RANRFAO (1 O.SYST, JK.L) 

call  RAMREAOI 1 1 . INF.MAXLEN.L) 

323  CONTINUE 
ipRsTPR.l 

IF(ACTUA)  (L) )  1.1.14 
34  IF  dNF(M)-  LEN)  3?, 51, 51 
51  CONTINUE 

1F(SYST(M,LEN))  1?.3?.33 
13  1hOLD=0 

lF(lr)ISTON)41  .41  .60 

C  OTSTRO  maintains  MINIMUM  LEVEL  IN  TqUR  AREAS* 

40  IFINCRNOOF (Kl) ,EQ.1  ,AN0.  LEVEL. EO.1)  60  TO  61 
iF(MINd)  .EO.li  GO  TO  6? 

CALL  MINIMUM! I ,K1 ) 

61  IF (SYST(M,LEN) )4?,6?.41 
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SYST(M,|  F  M)  =SYST  ♦iHOLf^ 

fin  Tn 

^,3  sonrM/i?.n»niiT  (i.) 

TF(SOO-n  7.7,H 
8  son=i,o 

7  lFr  =  «;YST  (M.LFN)  *SOO 

<;yst<m,lfn),=syst(m,lfn)-tff 

TTTxIFF^TT’ 

C - WHFM  PERCFNTAGF  was  applied  to  )  FNtTTY— THE  ENTITY  WAS  LOST  TO  SYST 

KOMaSYST (M.lFN) 

IF  (SYST(m,LFN)  .FO.l)  00  TO 
TFTPEPOMTdP)  .OT.O)  on  TO  1? 

PPROKT ( IF  I =1 ,00 

12  SyST<M,|FN)  =  SYST  (M,l  FN)*Pfc'ROlJT  (TF» 

C  PFTaTN  personnel  NEFOEO  for  MTNTMIIm  fill 
26  KOMsKOM-SyST <M,lFN) ♦ IHOLO 
Tr (SYST (M,LEm) -TOUT ) 3,4«4 
1  Ir.AlN=lr,ATN*SYST  (M.LFN) 

TOllT=IOltT-SYST  (M,LFN) 

TS=SYST (M»l FN) 

r  TF  KOM.GT,  O.PFRSONNFI  RFrAIN  in  SYsT(M,LFN>  AN5  INE(LEN)  RpmaTNS  TwF  Same 
IF(K0M,GT.  0)  00  TO  64 
TNF (LFN)=INr(LFM)-l 
64  TAI  I  RANWPTTF ( 1 1 , INE«VAXLFN,L) 

SYST (M,l 6N) =0 

«:YST  (M,LFN)  =SyST  (M,LEN)  ♦  <0M 
TALL  RANwRITE ( 1 n.SYSTi JK,L) 

OKKaM* 1 

fall  RANRFAO  { 10»SYST. JK«ICC) 

CALL  RAMPFAPdl  ,INF«maxLFN,ICC) 
r 

r - TyPF  3  INTO  (1,1) 

LEN2=1 

GOTO  (1100,1100,1300.1^00)  ITYPE  (TEND) 

1300  SyST(LEN,1)  s  SYSTILEN,!)  ♦  IS 
T.l.l*  LFN 
I  FN?s1 

IF  (TNE(LFN)-LEN2)  1  354 , 1  3S5, 1 355 

1354  TmE(LFN)=1 

CALI  RANWRITEdl  ,INE,maXLEN,ICC) 

1355  Gn  TO  25 

r 

C - Type  4  INTO  (1,1) 

1400  SyST(1,1)  =  SYST(1,1)  ♦  IS 
LEM?=TJJ=1 

Tf  IlNEdlrLEN?)  1454.145S.14S5 

1454  lNF(l)a  1 

Call  ranwritf  o i.ine.maxlen,  ico 

1455  GO  TO  25 
C 

C - TyPF  1  AND  2  (1  ,T) 

1100  IF  (M-MAXLFN)  24.24,23 

23  SYSTd.M  )»SYST(1»M  ).IS 
T.JJ'I 

LEN2sM 

iF(lNEd)-M)  S2,S2,b1 
5?  TMFd)=M 

CAt  L  RANWn  r  TF  d  1  ,  INF,MAXI.EN,1CC) 

53  continue 
-  60  TO  25 

24  SyST(1,(En)=  SYST  ( 1  ,LEM)  *15 
IJJ=1 
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;*>  o 


1 


irMPsL^N 

lF(Ir^F(n-LFM2> 

64  TMF(n«|.FN 

C«LL  RANWPTTF(n,INE.^'AXLEN»ICC) 

5S  CONTIMUF 
25  CnwTlNUF 
TP=TP»1 

M  pRTi^T  ■  M-1  maJ  Only 

UFM.S 

PRINT  loni,  K1,K3»I,J,I*;.IFF.IJj.LEn2«IP*ITeM 
Call  pack  (ki  ,ki.i»j,t<;,iff»  ijj,lEn2.ni »m2.n3  > 

CALL  RANwRlTF(nfNl,:?,IP) 

CALL  RANWRITE  ( 1 0 ♦SYCT , JK « ICC) 

r 

3?  IF(LFN-I)  <),q,6 
9  If(M-MONTH>  1,1,10 
10  M=M-1 

RO  TO  1?I 
f)  LFNeLPN-1 
60  TO  2 
C 

4  l6AlNsl6AIN*I0l)T 

CALL  RANWRITE(  1 0 « SYST , JK »L ) 

TALL  PANPEAO  I1C«SYST,JK.ICC) 

Call  RANRFAfK  1  I, 1NF«MAXLEN»ICC) 

r 

ROTO  (11  10»1  110,1310.1410)  ITYPEdFMO) 

C - type  3  INTO  (1,1) 

1310  SyST(LFN,1)  =  SYST(LEN,1)  ♦  lOUT 

-  IJJ«  LEN 

LEMZsl 

IF  (INE(LEN)-LEN2)  1357,22,22 
13R7  1NF(LEN)*1 
GO  TO  22 

- -Type  4  INTO  (1,1) 

1410  SyST(1,1)  »  SYSTd.D*  TOUT 
LEN2=IJJ»1 

IF  (TNEd  )-LeN2)  1457,22,22 
1457  INE(1)»  1 
RO  TO  22 

C -  .  . . . 

c - type  1  AND  ?  (1,1) 

1110  IF  (S-MAXLEM  21,21,20 
-  21  SYST(l,LEN)«SYST(l,LEN)4lOUT 

IJJ=1 
LEN2sLEN 

- lFdNE(l->-tEN)  57,22*22  -  - 

57  ImE(1)tLEN 
GO  TC  22 

- 20  SYST(1,M  )*SYST(1,M  )*I0UT 

IJJ*1 

LEN2*M 

- TF4lNE(I)-M)  5(S,5ft,22  - -  -  - 

5R  InE(1)sM*1 
22  CONTINUE 

-  IR»IP^1 

MpPInT«M-1 

ITEM»4 

- PRINT  loot*  Kl,K3,I,J,IS,IFF,iJJ,LEN2,  l»»^,  item  - 

CALL  PACK  (K1,K3,I,J,  lour,  IFF.UJ.i  EN2,N1  ,n2,N3  ) 

CALL  RANWRITF(13,N1,3,iP) 
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r»  r> 


I 


Call  RANWPlTF(n,lNC«MAXLEN«ICC) 
TALL  RANWRItE  (1 0*SYST, JK. ICC) 

- TALL  RANREAn  (lO.SYST, JK*L) 

SYST(M,LEN)=KOM-IOUT4SyST(M,LEN) 
call  RANRFAD( 1 1*INE,MAXLEN*L) 

-  Tr(lNE(LEN)-M)  58»S9*5Q 

*>8  ImE(LEN)xM 

Call  ranwpite(ii*ine*maxlen*icc) 

- RaNWRITE  (10*SYST.JK.L) 

OKK=M*l 

TOUT*0 

-  1  TFFSITT 

IF  (TfiRAOE(lENO)  .LE.O  )  fiO  TO  15 
PRINT  lon7«  TGRADEdEND)  ♦  lEND 

-  -  Call  TRN5P0SE 

15  CONTINUE 


LOCK  FORMATS 
lOni  format  <6H  lOCK)  lOTlO) 

100?  format  (6H  L^>CK?  10110) 

1003  FORMATd  IHOiniSTON  =  .12) 

1005  Format  (sh  tgraoes  ta.  6H  ifno=  U) 

1007  format  (SH  IGRAOEs  T4.  6H  IEN0=  U) 

return 

END 


FORTRAN  OTAGNOSTIC  RESULTS  poR 


NULL  STATEMENT  NUMBERS 
1003  1002 


LOCK 
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31/3?/330o  fortran  l3.1)/ySOS 


-  SUBROUTINE 

ElLUP 

COMMON 

NOEIRST 

, NOLAST 

•  NOTT 

.ITRtIME  ,INTEG 

COMMON 

IP 

,InE(48) 

,ICMD 

,OUTtT(?oO) 

- COMMON 

OUTST(?0ft) 

,OuTTO(2ftO)  - 

iWTjOIpno) 

COMMON 

PCT (?00> 

,PERO(?00) 

» Jif 

,SYSt (4R.48) 

COMMON 

PAdOO) 

,NAVd00) 

,LrNGTH(100)  ,OUT(lnO) 

COMMON 

INTOUR (??0 

) 

,ImSU8(220i  ,OUTTOUR(220) 

COMMON 

k0n(2?0) 

,Ol|TSUR(??0» 

,AFTcR(2?0) 

COMMON 

PER(220) 

,RATEdOO) 

.NEEDSdn)  «NE(10.10) 

COMMON 

lOSdO) 

,NFndO) 

iNFEdOOT 

,irtL;L  jUAXLEn 

COMMON 

NTOIIR 

,NP 

,ClOS 

•PERoTO(70rt) 

common 

POW 

,LSTRSTOdO) 

.LSTftSTTdO) 

COMMON 

LSYRSTSdO) 

,ITT 

,«axsur  ,  act  do) 

COMMON 

NSITdOl 

,lNOSdA,10) 

,1GRaOE(220) 

COMMON 

REP (220) 

,ITYPE(220) 

,PERoUT(220) 

-  —  -COMMON 

lOISTON 

♦  iS'iM 

,ORP5UMdOO»  ,ERlOdftO) 

COMMON 

HE6R0W (100) 

,En0R0w(100)  ,eEC,COL  (100) 

COMMON 

ENOCOL(IOO) 

,MaTSUM(100)  ,MAT6RDSlinO) 

COMMON 

Type  doo) , 

SUB (100) 

,NuM(100) 

, ACTUAL (1  on) 

COMMON 

nprlev 

,nt 

,IhOL0 

,LEN  , LEVEL  tM 

COMMON 

NCPNOOFdOO) 

,MtN(100) 

,6RPtNPRrlOO) 

COMMON 

MAXOEPLdOO) 

Integer 

Type 

,Surt 

, Actual 

,ENOfOL 

integer 

CIOS 

,SYST 

,OuTSUR 

,OUTtOUR  , after 

Integer 

ACT 

,OUTTO 

,0uTS0 

,OUTfT  ,OUTST 

InTFGER 

PERU 

,GRPSUM 

,BFGRO)i( 

,ENO«OW  ,BE6C0u 

Integer 

PEROTO 

INTEGER  OK, 

OKK 

Ml*N?*N3*0 
TF(tFIl.L-l)  S,in.lO 
5  on  15  I*1»ITT 
IS  NF(I)«0 
10  lrOUNT*0 

— —ChFCK  each  tour  to  FTIL  PFOUIRFMF.NTs  ACCOROTNG  tO  FILL  HllLFs, 
DO  ?0  TFnO=1«NP 
Kl.TNTOURdEND) 

K?»OUTTOUR( IFNO) 

K3=TNSUB( IEmD) 

KA«OUTSUB( IENO) 
nK=o 

TF(lFlLt-n  35,40,40 
40  lFTNn*NF(Kl,K3) 

60  TO  4S 

35  lFlNn«NEEOS(Kl) 

45  irr«(Ki-n*MAxsuR»K3 
IF (PFR{ TFMO) )  75,P6,30 
30  iFlNOxPFRdFNO)  *  IFIM) 

?f.  CONTINUF 

TFdFlNO)  140,141,14? 

140  iFlNOeO 

141  60  TO  ?0 
14?  CONTINUF 

iFlOlJTTOURdFNO)  »  R0,Pn,7F 
75  !r«(K?-l)  *  MAX5MB^K4 
C  INITIALLY,  SYSTEM  INPUTS  THOSE  WHO  HAVE  COhPLFtFO  THF  TOUPs, 

LA*LFN6THdC)  ♦! 

CALL  RANPFAf)dn,SYST,,)K,  (Cl 
no  ?5S  TMsI ,MAXLFN 
IF(5YST((  A,IM) ,LF,0)  00  TO  ?55 
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TF (prROiiT (  IEmD)  ,r,T,0)  Rn  to  5S 

pfrOijt  (  *1  .on 

•yty  K0M««;Y5T  «LA,1M)  •PFROIIT  ( rFKO» 

IF  (SYRTd  A.  IM)  .FO,  1 )  KOM*! 
lF(t<nM)  ?A,?,?A,?,25h 
256  If  ( TFlNn-KOM)  fo,60»A>F 
fcO  rONTIAJllE 

— — ni'mbfr  availahlf  is  not  Sufficient  to  fill  rf.outrfmcnts, 

SyST(|  A*IM)>(;ysT  (LA«TM)-IFIN0 
CALL  RANPEAn( lO.SYST, JK*ICC» 

POINT  UT.SYST  (LA.IM)  ,TM,LENr,TH(lC)  .KOM 
S-RTCUTN  )»SYST(1.IM  )*IFtNn 
IP=IP*1 

PRINT  InnS.Kl ♦K1.K?,K4«IFIN0»Im,TP 
CALL  PACK  (KI,Kl,K2tK4«lFlNn,0«TM,l«|yi1  .N2«H3t 
CALL  RANwRiTEt  I3«NI«3«lP) 

CALL  RANWRTTF(10*SyST,JK«ICCI 
OKslFlNn 
iFlNOxO 

ArT(Ki)*ArT(Kn*oK 
actual ( ICC) sATTMAL ( ICC) ♦OK 
no  TO  70 
65  CnwTTNUF 

- UjjMBFp  available  ExCFFn<;  rEQuIrFMENts. 

TFlNn«TFTNO-KOM 
OKtKOM 

SyST  (LA,  IM)  *«;y«;T  (LA,  1H)-K0M 
CALL  RAWRFAn(IO,SYST,jK,ICC) 

SVST(1,TM  )sSYST(I,lR  )40K 
TP*IP4l 

PRINT  lr,n4,Kl  ,Kl,K?,KA,nK,IM,IP 
CALL  PACK(K1 ,K3,K?,KA,OK,0«IM,1,NI,N2»N3) 

CALL  RANwRItE (13,Nl,3,IP) 

CALL  RANWRITE(10«SYST,JK,ICC) 
act (K1 ) *ACT (Kl) ♦OK 
actual ( ICC) sACTUAL ( ICC) ♦OK 
0KKsM4 I 
00  TO  250 
255  continue 
258  CONTINUE 

TF(  IFTNO)  2f,0,70,?F.? 

260  IfINO=0 
00  TO  70 

C  THERE  SHOULO  NOT  0E  ANYONE  BEYOND  THE  TOUR  LIMITS  OR  LENGTh(IC). 

262  IF (AFTER! lENO) -length ( IC) >50, 50, 70 
50  I0UT*TETND 

TL0SS«TGAIN*IEE30 

MMMMsMAXSUR 

“MYaHiXStlR 

IF«IEN0 

C  CFARCH  EOR  personnel  IE  THEY  ARE  AVAILABLE  aBOVe  MINIMUM  LEVELS,  FOR 
C  CRTTtCAL  TOUR  LEVELS  IGNORE  MINImuM  LEVEi  S. 

TFdniSTON.EO.O)  GO  TO  ?S 
IF (MIN(K?) .EO.O)  GO  TO  ?5 

teilevel  .Eo.i  .and.  ncrnodE(ki»  .EQ.i)  go  fO  ?5 

GO  TO  20 

25  CALL  L0CK(K2,K4,AFTER(IENn) flLOSSflBAlMtlOUTflFr.My.ICCtlE.KltKai 
Print  loo?,  K2,K4,AETFR(lENd»,fL(5SS,T(iAIN,IOUT,TFE,|iY*ICC,TE,*<I,K3 

MAXSURaMMMM 

aCT(K2)sACT(KP)-IGAIN-IFE 
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OK=TRATN 

TFINDsIOIIT 

- ^kCTUALi  ir>*ACTOAL(IC)-IGAlN-IFf 

act (Kl) =ACT (K1 ) ♦OK 
actual ( ICC) = actual (ICC)  ♦OK 

-  .|f.<^FH_L-l)  !^0.151*I51 

ISl  NF (K?,K4) =MF (K?,K4)^IGATN4!FF 
ISO  needs (K?) =NEEOS(K2) ♦TQAIN^IFF 
-  -70  TF(IFILL-I)  fl5»<»0»90 

90  NE(K1,Ks)=IFIN0 
P5  NEEDSIKI )=IFTNn 

-  -  00  TO  i?n 

C 

C- - IMPI'T  SOUGHT  From  outside  SYSTEM, 

-  «o  Cali  outstof.  ( ifno) 


V  0 
-  ?0 
c 
c 

1001 

100? 

1005 

100#> 


00  TO  l?0 
IcOUNTsTCOUNT^OK 

continue 

FTLUP  formats 
F0RMAT(SH  FIlUP,4I10) 
format  (l?H  BEFORE  LOCK  1218) 
format  (12H  AFTER  LOCK  12IR) 

Format  (ph  fIi.upI  sir,8x*2T8) 
format  (8H  FiLUP?  SIR, 8X, 218) 

pfturn 


fmo 


FORTRAN  OTAGNOSTIC  RESULTS  pOR 


NULL  statement  numbers 
1001 


FILUP 
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ri  n 


3l/3?/330n  FORTRAN  (3.1)/mSOS 


SUBROUTINE 

TIMFIIP 

.  - - 

- — -  —  - -  - - -  .  -  .  . 

COMMON 

NOFTRST 

«NOLAST 

♦NOTT  ♦ITRtIME  tlNTEG 

COMMON 

IP 

,INE(48) 

♦IFNO  .OUTtT(2oO) 

-  -  f^MMON 

OIJTST  (200) 

.OijTTo(?ob)  '  .OuTsOCZnO) 

COMMON 

PCT (200) 

,PFRn(200) 

♦ Jy  .SYSt (48>48) 

COMMON 

Pa ( 1 00) 

»NaV(V00) 

♦Length! ion)  ,out()oo) 

-  ■ • COMMON 

INTOIJR  (220) 

♦  InSUB(220)  tOUtTOriRC^Pn) 

COMMON 

100(220) 

,OljTSUR(220)  .AFTpR(220) 

COMMON 

PER(220) 

»RaTE(100) 

.NEEpS(lo)fNt(10,10) 

-  ■  COMMON 

lose  10) 

♦NFOCfO) 

♦  Ne'e  (1 00)  ■.  1E!lL  tMAXLEN 

COMMON 

NTOUR 

.NP 

*ClOS  ♦PERdTOI20o) 

COMMON 

P(5W 

♦LSTRSTO(IO)  .LSToSTTdO) 

- COMMON 

LSTRSTS(IO)  ^ 

♦ItT  ♦MAXsUB  tACTClO) 

COMMON 

NSIT ( 10) 

•lNOS(lO,10)  .I6RaDE(22o) 

COMMON 

REP(220) 

.ITYPE(220)  .PERoUT(220) 

- — n^-OMMON 

'TOISTON 

dSUM 

♦6RPSUM(i0n)  tPRlOCioO) 

COMMON 

RFGPOW (100) 

♦EnDROW(IOo)  »REGCOl(100) 

COMMON 

ENOCOL (100) 

♦  M(»TSUM(10n)  _  jMATGRpSdOO) 

-  rOMMON 

TYPF(IOO) « 

SUR(IOO) 

♦  NiiMdOO)  ♦ACTgALdOn) 

COMMON 

nprlev 

♦  NT 

♦IwOLD  ♦LEN  fLEVEL  tM 

common 

NCRNOOE( 100) 

♦  MtN(IOO)  ♦GRPTNPRdnO) 

—  --  COMMON 

MAXOEPLdOO) 

TmTFRFR 

Type 

♦  SUB 

♦Actual  ♦EnocOL 

Integer 

CIOS 

♦  SYST 

♦OuTSUB  .OUTyOUR  ♦AFTER 

-  Integer 

ACT 

♦  OljTTO 

♦OuTSO  ♦OUTtT  »outst 

Integer 

PFRO 

♦GRPSUM 

♦begrow  ♦ENDrOW  ♦BEGCOl 

Integer 

PFROTO 

'  OTMENSTON  SYS(4«»Afl) 

TMTFfiFR  SYS 
PRINT  n«TTRTlME 

jk*haxlfn*maxlem 

TFdTRTIME.EO.?)  GO  TO  10 
TSTARTs] 

TSTOPsNOFTRST 
on  To  n 

IP  T«;TAPT=N0FIRST*1 
TSTOPrNOTT 

1  1  00  1  TsTSTARTdSTOP 
NlsN?sN3*n 

K 1  sOlJTTT  (  T  )  SK2=0mTT0  (  T  )  fK^rOUTST  « I )  *K4=0UTS0  ( 1 ) 
TF(nHTTT(T)-l) 

transfers  RFHAvF  LIKF  LOSSES  TO  SyStEM, 
f,  K=  (ntiTTO( n -1 )  *  maxsijr*outso(I) 
call  PANRFAOC 10,SYST.JK,K) 

MSPaPF.ROC  I ) 

OO  R  L=1 .MAxLFN 
TF(SYST(MSP,L) )  8f8,R 
q  NMN=SYST(mSP,L)»PCT(I) 

SYST(MSP«L)=SYST (MSP»L»-NNN 

PRINT  3.  nUTTT (T) .OUTS f « I ) .OUTTO Cl ) .OUTSOCI) .PCt (I)  ,PERD(I) , 
VSyST(MSP.I  ) 

IP=IP*1 

CALL  PACK (K1 ,K3,K2,K4.0.NNNfMSP*L  *N1.N2.N3) 

Call  PANWRTTE(n»NJ,3.TP) 

8  CONTINUE 
jRO  TO  1 

C - TRANSFERS  ARE  FROM  ONE  TOuR  TO  ANOTHER  wITHTN  TWE  SYSTEM, 

7  JsCOlITTT  ( I)-l)  *MAXSUB»0UTST(I) 
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Ka(OIITTn(I)-l)*MAXS(jn*OUTSO(I) 

CALL  RAMRfAn{  in,SYST,.JK»K) 
tall  PANRFAn< lO.SVS, JK,J) 
mSP=PERI)(I) 

DO  ?  L=I.MAXLFN 

-  -  ff  (SYST(MSP«L))  ?»?»5 

s  cnMiiMUP 

NmN=SYST (MSP,L) *PCT ( T ) 

•^YSKMSP.DaSYST  (MSP,L)-NNN 

MSPTOsMSP 

lF(PERnTO<T) ,GT.O)  MSPTn=PFRDTO ( I ) 

-  -  sys<mspto,l)=sy‘;(mspto,l) ♦nnm 

PRINT  nilTTT(T),OUT«;T{n*OuTTO(I).ni)TSO(H.PCT(I).PERO<I)  , 
ISVST(MSP.L) »SYS(mSPTO,H 

|palf>*\ 

CALL  PACK (K1 ,K1,k2,K4»MNN,0*MSPT0*L.M  »N?«N1) 

Call  PANWRlTF(n,Nl,1,TP) 

-  ?  continue 

TALL  RANwRITF (10,SYS.JK«J) 

1  CALL  RANWRlTE(ir,SYST.jK,K) 

TTMFiJP  FORMATS 
1  format (4I  1  n,FH. 4.41  in) 

11  FORMAT(7HnTlMEl)P.?X,  lOHITRTlMF  =  .1/,) 

return 

ENO 

FORTRAN  OIARNOSTIC  RESULTS  FOR  TlMEUp 


NO  ERRORS 
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31/1?/330o  FORTRAN  (3,1)/mSOS 


COMMON 

COMMON 

'TTVWMON 

COMMON 

COMMON 

“COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

"COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

ImTFGFR 

InTFRFR 

TNTFOFR 

intfgfr 


OUTSTDE(I)  -  -  - 

NOFTRST  »NOLART  *NOTT  ,ITRtIMF  *INTFG 

IP  »INE(48)  .IgMO  .OUTrKanO) 

OllTST  (?00)  «OuTTO(?rtO)  «OUl^ ^0(200) 

PCT(200)  »RFPO(?0ft) •Jk 

PA(IOO)  ,NAV(iOO)  tLENGTHdOn)  .OUTClnO) 

INTOUR(2?0)  «InSUB(2?0)  *OUTTOiiR(2?6> 


100(220)  ,OUTSUB(2?0» 
PFR(220)  ,RaTE(100> 

105(10)  tNFO(lO)  ♦Nef(IOO) 
NTOiJR  «NP  ,ClOS 

POW  .LSTRSTO(IO) 

LSTRSTS(IO)  ,Itt 

N5IT(10)  «INOS{10*10) 

REP(220)  «ITYPE(?20) 
lOTSTON  dSUM 
RFGROW (100) 

ENOCOLdOO) 

TYPF(IOO)  .RUPdOO) 

NPRLEV  ,nt 

NCRNOOEdOO) 

MAXOEPLdOO) 


,AFTeR(220) 

.NEEpiSdO)  *NEd0,10) 

•  IFI^L  *MAXLEn 

,PERqTO(200) 

.LSTpiSTTdO) 

*MAX5UB  fACtdo) 

.I6RaOF(226) 
,PERoUT(220) 

•6pp5UM(]0o)  •PRTO(ToO) 

♦  ENDROWdOo)  •BEGCOLdOO) 

♦  MATSUMdOO)  tMATGRpRdOO) 

♦  NuMdOO)  .ACTuAL(TOn) 

.LEN  .LEVEL  .M 

.GRPrNPRdoO) 


•Ihold 
•MtN(IOO) 


Type 

CTOS 

ACT 

PFRO 

PFROTO 


•  Sim 
.SyST 
.OUTTO 
,GRPSUM 


T^'TFGFR 
TNTFGER  OK 
Ml=N?sN3=0 
KAsOOTSUBd) 

K3s0UTT0(IR(  T) 

K2=iNSimd) 

Kl=TNTOUP(I) 

C  IFkAsO  this  means  that  PROGRAM 

C  T0S(K4)  a  VECTOR  OF  NEW  INPUT 

IF(TFlLL-n  S,SdO 
IfTNo=nE(K1 ,K?) 

GO  TO  P7 
iFlNOaNEEDSlKl ) 

If ( Tftno.lE.o)  go  to  Rl 

TF(K4-1)  R2,R3,q3 
TF(RFf’d)-l.O)  84»fl5«05 
SSAVFalFiMf) 

TFlNOaPnw*SSAvE 

GO  TO  PA 
iFiNOrSSAVE 
KKalOSd) 

TOSd  )  =  IF  INQ 
KjaK4 


♦Actual 

♦OyTSUH 

♦OUTSO 

♦begrow 


.ENO^OL 

.OUTtOUR 

.OUTtT 

.ENOrOW 


♦after 

•OUTST 

♦BEGCOl 


IS  USING  VaRIArLE  INPUT 
during  this  timc  period 


10 

s 

07 

8? 

04 


RS 

Rh 


K4a1 

83  IF(PFR(T))  20. 20. HI 
81  OKr lOS (K4) *PEP d ) 

Ir(lOS(K4)-  TFTNO)  02.31,3) 
92  iFlNOalFINO-OK 

IOS(K4)rTOS(K4)-OK 
GO  TO  3S 

20  If  dOS(K4)-lFlNO)  25,30.30 
25  lFlNn=IFlN!)-10S(K4) 

Ok=  T05 (K4 ) 

TOS (K4) =0 
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o  n 


r,0  TO 

30  Ok-IFTNO 

- 31-In5;<K4)  ■  I0S<K4)-IFIND 

nK=iFiNn 

TFInO=0 

-  . . {F(IFTLL-I)  40,«i0*50 

*=.0  NE(K1  ,K2)  =NF.  (K1  ,K?)-OK 
40  NEEDS(K1)=NFEDS(K1)-Ok 

- -  TN0S(K1  ♦K?) =Ok 

ICC=(K1-1  )*MAXSIje*K2 
CALL  RANREAn(10.SYST.JK*ICC) 

-  SYST<l*n=SYST(Ul)40K 

CALL  RANWRITE(10»SYST#JK,ICC) 

TP=TPH 

-  - PRINT  loos*  K1  *K2*K3*K4*0K*  IP 

call  Pack (Ki *K?*K3*K4.0K,0*1 *1,N1*m?,N3) 
CALL  RANWRITF.  ( I  3*N1 ,3*  TP) 

-  ACT(Kn*ACT<K|  )  *OK 

actual (ICC)=ACTUAL(ICC)*0K 
TFiKJ-n  R'',8Q,88 
8R  SSAVF-IFIMP 
TOS(n=KK 
KA  =  KJ 

0B  CONTINUE 

CALL  RA^JPEAO(1  1  *rNE*MAXLEN.ICC) 
TF(TNF(1)-1)  R0.91.QI 
VO  TmE(1)=1 

TALL  RANWMITE(11.TNE*MAXLEN.ICC) 

9l  continue 

OUTSIDE  formats 

1005  format  (12H  OUTSIDEr—  1?IB) 

return 
EmO 

FORTRAN  OTAGNOSTTC  RESULTS  POP 


NO  ERRORS 


OuTSIdE 
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3l/l?/330o  FORTRAm  (3.1)/mSOS 

- SUBRn'ITTNE  TRNSPOSE 


COMMON 

NOFIRST 

.NOLAST 

.NOTT 

.ITRtIMf  .integ 

COMMON 

IP 

*INE(48) 

«  ICNO 

.OUTfT (29O) 

COMMON 

OIJTST  (200) 

«OuTTo(26o) 

.OUT50(2p6> 

COMMON 

F'CT(200) 

♦PFRD(200> 

»  Jk 

.SYSt (40. AS) 

COMMON 

PA(IOO) 

.NAVdOO) 

•  LENGTHdOO)  .OUTdoO) 

COMMON 

INTOIIR  (220) 

♦lNSUR(220)  «OUtTOuR (?20) 

COMMON 

100(220) 

.OuTSUR(220> 

, AFT^R (220) 

COMMON 

PER(220) 

.RATEdOO) 

.NEEnS(lo)  .NE  do.  1  0) 

COMMON 

IOS( 10) 

.NFDdO) 

»NeE(100) 

.IFIlL  .MAXLEN 

COMMON 

NTOUR 

,NP 

•  ClOS 

.PERpTO (200) 

COMMON 

POW 

»LSTRST0{ 10) 

.LSTpSTTdO) 

common 

LSTRSTSdO) 

♦  ITT 

.MAX5l)R  .ACfdo) 

COMMON 

NSIT ( 1 0) 

.INOSdo,  10) 

.  IGRaOE (226) 

COMMON 

RFP(220) 

,lTYPc(?20) 

.PER(5UT(920) 

COMMON 

lOISTON 

♦  ISUM 

.grpsum(io6)  .l»Rio(vo6) 

COMMON 

HEGROW ( ino) 

.EnOROK(IOo)  *BEGCOl(1()0) 

COMMON 

ENOCOL (100) 

.MATSUMdOO)  tMATGRPS  (Too) 

COMMON 

TYPE(1oO)» 

SUBdOO) 

♦  NuMdOO) 

.actual (lo6> 

COMMON 

nprlev 

,nt 

♦imold 

.LEN  tLEVEL  .M 

COMMON 

NCRNOOEdOO) 

♦  MTNdOO) 

.GRPiNPRtloO) 

COMMON 

MAXOEPLdOO) 

IMTFOFR 

Type 

.SUB 

♦Actual 

.ENOcOL 

I nTF  GPR 

CIOS 

.SyST 

.OuTSUR 

.outtour  .after 

intfgfr 

ACT 

.OUTTO 

♦OuTSO 

.OUTtT  *outst 

InTFGFM 

PEHO 

.GRPSUM 

.begrow 

.ENOrOW  .BEGCOl 

Intfgfr 

PEMOTO 

ir  (IGRADEdENO)  .LF.  0 

)  GO  TO  R99 

ENTRY  FLIP 

no  10  lai. 

MAXLEN 

on  in  Jal, 

MAXLEN 

ISAVE*  SVSTdfJ) 

SyST ( I,J)r 

SyST  (.).  I  ) 

10  SyST ( J«I! =  ISAVF 
ReTURM 
END 

FORTRAN  DIAGNOSTIC  RESULTS  FOR  TrNSPoSE 


NO  ERRORS 
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31/l?/130o  fortran  O.D/mSOS 

SliBRODTTNF  SiiMMARY(NIIMElFm,  INf)IV«INTnT) 

C  APRIL  lQ<^R-PROr,RAMMER  WllT 

C  SI'MMARY  SltBROUTlNE-OHTATMS  SUMS  OF  PARTIAL  aND/qR  COMPLETE  mATPiCE 

r  S  AMO  OROllPS  OF  MATRICES 

C  SIIBROUTINF  parameters 

r  IF  TNDTV  .OT.  0. TOTALS  FOP  WHOLE  MATRICES  ARE  qRTATNED.  VALUE  OF  tNoIw 

c  oeterminfs  matrix  to  be  summeo 

r  tntot  contains  matrix  total 

r 

COMMON 

common 

COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
InTFGFR 
InTFGFR 
InTFGFR 
InTFGFR 
INTFGFR 
C 
C 

c  Section  \ 

C  SUMS  INOIv/IOilAl.  matrices 

c  NuMELEMsNmmHFM  OF  elements  IN  MATGRdS  VECTOR 

c  maTGRPSsVFCTOR  which  oeterminfs  which  matrices  oR  vectors  adf  TO 

C  HE  SUMMED  and  the  aggregate  SUMS  TO  ME  O0TAInEO 

r 

IF (IndIv.GF.I )G0  to  4 
print  inn 
PRINT  ino 
4  istartsi sisuM=n 
j=n 

no  in  t  =  i.numEi,fm 

MATS||M(I)  =  n 
iFdNOlV.GF.l  )GO  TO  ? 

IF  (MaTGPRS  (  n  .FO.OGO  TO  S 
J=J*1 

maT=matgrps ( T ) 

IBFGR=HFGR0W( J) 

IENDRsFNOPOW ( J, 

TBEGc=eEGrnL { J) 

TENDCxEndcOI.  ( J» 

GO  TO  1 

?  I0FGR=IBEOC=1 


NOFIRST  .NOLAST  ,NoTT  .ITRtIMF  .INTEG 

IP  ,1NE(4B)  «lENr)  .OUTiTiznO) 

OUTST<200) ,OUTTO<?00»  .OUTcOlZOO) 

PCTIpno)  ,PF«O<?0n»  *Jw  .SYSt  (aS'.AB) 

pa(ioo)  tNAVdoo)  •LcNGTH(inn)  «OuT(inn) 

lN70UR(2?n)  .lNSUR(22n»  tOUTTOuR ( 2?0 ) 

IOO(2?n»  ,OuTSUB(2201  .AFTfPIgaO) 

PFR(2?0»  ,«ATE(10n)  .NEEpS(lo)  tNEdn.lO) 

lOSdO)  .NFOdft)  »NEEdOO)  ,1E1lL  aMAXLEn 

NTOUP  ,NP  «ClOS  .PERdTO(200) 

POW  .LSTPSTOdO)  .LSTrSTT(10> 

LSTRSTSdOl  aItT  .MAXsUB  tACTdn) 

NSIT(lO)  .iNOSdOtlO)  .I6RA0E(22ft) 

RFP(22ni  .ITYPEC2201  .PERoUTl??0) 

lOTSTON  ,ISUM  ♦GoPSUMdOO)  •PHfOdOO) 

HEGROWdnnt  aEnoROWIIAo)  •BEGCOLdOO) 

ENOCOLdOnj  .MATSUMdOO)  •MATGRPSdOO) 

TYPEdAnt«StJB(lOO)  tNuMdOA)  .ACTuALdAn) 

NPRLEV  »NT  aIwOLH  ,LEN  aLEVEL  »M 

NCRNOOEdOA)  aMtaKIOA)  .GRPrNPRdOOl 

MAXDEPLdAO) 

Type  ,suh  «actual  .endcol 

CTOS  .SyST  aOuTSUB  .OUTtOUR  aAFTER 

act  ,0UTT0  .OuTSO  ,0UTtT  aOUTST 

PERO  .GRPSUM  aHeGROW  .ENOrOW  aBEGCOl 

PERDTO 
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or>  !->  r>o">onr>r» 


lENnRsLFNfiTHdNllIV) 

IENDCsMaKLEN 

MiTsTNDTV 

1  tall  PANPFAOdO.SYST.JK.MAT) 
no  1  Kr THFhP. lENOR 
no  1  I  =lMFGC.IENnC 
I  MATSljM(  T  )  =MATSllMd)  ♦SY*;?  (K«L) 

If (iNOiv.GE.i )r,n  to  \\ 

PPINT  10 1. INTER, MATGRPS (1 ) »MATSUM(1> ,lHEGRtTENDRtlBEGC»lENDC 
60  TO  10 


section  ? 

SUMS  GROUPS  OF  MATRICES  ANO/OR  VECTORS 

OrPSUM«SUmS  of  groups  of  matrices  Or  VECTORS  DESIGNATED  IN  mATGRPS 
VECTOR 

S  ITENDzI-1 
I(;llMslSiiM,  1 
RrPS|JM(ISUM)=0 
no  n  KalSTARTdlFNO 
8  RPPSlJM  (  T  SUM)  =GRPSUM  ( I  SUM)  ♦MATSUM  (K) 

PRINT  102*  GPPSUMdSUM) 

ISTARTsT 
10  continue 

CoLL  PANWPTTE(1*,GRPSUM*IS')M,INTFG) 

Go  TO  1? 

U  !fiTOT*MATSUMd) 

12  (ONTTNUF 

SURPOUTTNF  SUMMARY  FORMATS 

100  format  (/lIHOOUTPUT  FROM  SUBROUTINE  SUMMARY/) 

101  FORMAT(IS,T«,n2,*X,I2,lH-,I2,3X,I3,lH-,I2) 

10?  FORMAT(4TX,I10) 

103  format (7HOPFrIOD,3X,AMNODE,3X,10HNOdE  TOTAL. 3X,aHR0WS»3X,7hc0LUMNS 
13*,13HCLUSTER  TOTAL) 

return 

End 


FORTRAN  OIAGNOSTIC  RESULTS  FOR  SuMMArY 


NO  ERRORS. 


rt  r^n  i  o  r>  «*> 


31/32/3300  fortran  {3.1)/mS0S 


-  PROGRAM  FILL  - 

COMMON  A(4R«4R)*  NB(4R«4A) 

COMMON  !HT(2000) 

- RFAL  PFR»N  - 

intfgfr  a 

C - Mr  matrix  SI7F(SQUARE) 

—  .  RFAO  10]  «M  - - 

i«;tzf»m»m 

c - I«;TART  s  starting  MATRIX - ISTOP  *  ENDING  MATRf X ( IwCLUSIVEl 

- V  RFAD  ini»  ISTART,  ISTOP^  - 

IF  <ISTART-1)  RqR»2«2 
?  DO  f*0?  lAs  ISTaRT«  ISTOP 

PRINT  107  - 

call  RANRFAO  (10,A,ISIZE*1A) 

READ  101«K,IZER0 

- K«  NDMRFR  of  data  CAROS  - 

- IF  IZERO  EQ  0  MATRIX  WILL  BE  ZEROEo  THEN  OaTA  ENTERED 

- IF  IZERO  IS  I  OATA  WILL  0E  EnTERFD  INTO  EXISTING  MATRIX 

-  IF  (IZERO-l>  3«4«4  - 

3  DO  10  1=1. M 
no  10  J=1.M 

10  Afl.J)  »  ft  -  - 

4  n=N«Ir  JrPFRxO.O 
no  ?0  KA*1 .K 

•  read  100.  N.T.J.PER 
IF  (I.LF.O)  GO  TO  20 
IF  (J.LE.O)  60  TO  20 

Na  NUMBER  1 4  FORMAT  -  -  - - 

la  STARTING  ROW 
Jr  STARTING  COLUMN 

C- - PER  a  «  OF  N  IN  EACH  ELEMENT 

C— — IF  PER  s  ft.  loo  PERCENT  IS  ASSUMED 
IF  (PFR.EO.O)  15.16 

-  15  PER  al  .  - 

16  11=100./  (PER*10ft.) 

DO  20  Lal.II 

na*  n»  per 

IF  (I.GT.M)  ?3,24 


24 

IF  (J.GT, 

M)  ?3 

.25 

25 

IF  (IZERO 

-1)  26.21.21 

26 

IF  (Ad. 

J))  22.21 

21 

A'(I.J)  a 

NA 

IaI*I 

JaJ*  1 

GO  TO  20 

22 

print  103 

.  I  .  J 

L  =  II 

GO  TO  20 

23 

PRINT  104 

.  I  .  J 

Lall 

20 

CONTINUF 

30 

Call  ranwritf 

(Ift.A.ISIZE 

.lA) 

PRINT  106 

.lA 

Call  RANRFAO 

(lO.A.lSIZE. 

lA) 

DO  31  1= 

1  .M 

DO  31  J=I 

.M 

ISAVE=  A( 

T.J) 

A<I.J)a  A 

( J.T) 

31 

A(J.T)=  ISAVE 

I  aT  =  o 

-  TO  - 


jTi*n 

no  #>nn  iIs1»tsI7f«m 

JTI*  JIT*1 

ITA= 

on  600  J»  JUUA 
TF  <A (J) )  601 *600*601 
601  IaT«  IAI*1 

TbI  (TAI-l)a  JTI 
IPTCIAT)*  A(.I) 
lRHTAI-?)a  J-J1*1 
600  CnNTTNUF 

ir  (TAT)  60?,60?*603 

603  PrTNT  604*  (THT(J)«  Jal.lAl) 

604  format  (7(4X*?I4«I7)) 

60?  CONTTMIJF 

60  TO  1 

C— — Data  format 

TOO  format  (F6,o.?t?*F2,?) 

101  format  (?I4) 

103  format  (14H  OATA  for  ROW  T2*  5H  COL  !2*/  36H  NOt  ENTERED  PRIOR  0 
XATA  STORFO  TMERF) 

104  Format  (14H  data  for  row  I2*5H  col  Ig*/  34H  ENtEREO  but  data  fxc 
XFFOS  matrix  ) 

106  format  (1X*T3) 

107  format  (///nX,  3HC0L  ♦AH  ROW  ) 

9R9  FMO 

FORTRAN  OlAfiNO'iTTC  RESULTS  FOR  FILL 


NULL  STATFMFNT  NUMHFRS 

....  3Q 

ERROR  TYRF  1046  OETECTEO  AT  1  STATEMENT  BEYONO  STATEMENT  NO,  22 

-TME  ITUNNING  TNOFX  IN  A  00  MAY  HE  CHANGED  WITHIN  THE  LOOd, 

ERROR  TYPF  1046  DETECTED  AT  I  STATEMENT  BEYoNO  STATEMENT  NO.  23 

TMt  RUNNING  Index  in  a  do  may  BF  changed  WITHIN  the  LOOP, 
COMPASS*P*L«X. 
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3300  COMPftsS/M<;OS-COMP/»SS  V 


XS‘»P 


entry  PAcK»iinPACK 


oonno 

01077777 

01 

0 

77777 

0 

PACK 

UJP 

•  ft 

00001 

47100107 

47 

0 

P00107 

1 

STT 

XRt  ,  1 

oooo? 

14600000 

14 

1 

00000 

7 

ENA 

OR 

00003 

147O000O 

14 

1 

00000 

3 

ENO 

OR 

00004 

54100000 

54 

0 

POOOOO 

1 

LDl 

PAcK»1 

OOOO'i 

?1 500006 

?1 

1 

00006 

1 
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